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PREFACE 


This report is a sequel to an earlier report titled ,! Aeroelast ic 
Effects in Multirotor Vehicles with Application to a Hybrid Heavy 
Lift System, Part 1: Formulation of Equations of Motion”, (NASA 

CR-3822, August 1984). 

The research effort reported herein was carried out in the Mechanical, 
Aerospace and Nuclear Engineering Department at UCLA by Dr . C. Venkatesan 
and Professor P. Friedmann who served as the principal investigator. 

The authors want to take this opportunity to express their gratitude 
to the grant monitor Dr. H. Miura for his numerous constructive comments 
and suggestions, as well as for much of the numerical data used in Section 
4.2.1. 
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SUMMARY 


This report is a sequel to the earlier report titled M Aeroelastic Effects 
in Multi-Rotor Vehicles with Application to Hybrid Heavy Lift System, Part I: 
Formulation of Equations of Motion". The trim and stability equations are 
presented for a twin rotor system with a buoyant envelope and an under slung load 
attached to a flexible supporting structure. These equations are specialized 
for the case of hovering flight. The stability equations are written in multi- 
blade coordinates. The total number of degrees of freedom for hybrid heavy 
lift vehicle consisting of two four bladed rotors is 31. Hence the stability 
analysis yields a total of 62 eigenvalues corresponding to these 31 degrees of 
freedom. A careful parametric study is performed, and used subsequently to 
identify the various blade and vehicle modes. The eigenvalues are identified 
by relating them to the physical degrees of freedom present in the system. 

This identification is based on a parametric study in which the fundamental 
parameters governing the system are varied. The coupling between various 
blade modes and vehicle modes is identified. Finally, it is shown that the 
coupled rotor/vehicle stability analysis provides information on both the 
aeroelastic stability as well as complete vehicle dynamic stability in the 
longitudinal and lateral planes. Also presented, in this report, are the 
results of an analytical study aimed at predicting the aeromechanical stability 
of a single rotor helicopter in ground resonance. The theoretical results are 
found to be in good agreement with the experimental results available in the 
literature, thereby validating the analytical model for the dynamics of the 
coupled rotor/ support system. 
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1 . INTRODUCTION 


This report is a sequel to the previous report entitled "Aeroelastic Effects 

in Multi-Rotor Vehicles with Application to Hybrid Heavy Lift System, Part I: 

Formulation of Equations of Motion" [Ref. 1], in which the equations of motion 

governing the aeroelastic behavior of an approximate model representing an Hybrid 

Heavy Lift Airship (HHLA) (Fig. 1) were derived. The equations derived in Ref. 1 

were representative of a somewhat simplified model shown in Fig. 2. The model 

consists of two rotors, a buoyant envelope and an underslung load, attached to a 

* 

flexible supporting structure. The various degrees of freedom, considered in 
deriving the equations of motion, are flap, lag, torsion for each blade, rigid 
body translation and rotation of the complete vehicle and the degrees of freedom 
representing the normal modes of vibration of the flexible supporting structure. 

It is useful to review some of the more important assumptions used in deriving 
the equations of motion, namely: 

1) The rotor consists of three or more blades. 

2) The rotors are lightly loaded. 

3) The rotors are in uniform inflow. 

4) There is no aerodynamic interference between the rotor and the buoyant 
envelope. The aerodynamic model used for the rotor blade is the quasi- 
steady aerodynamic model with apparent mass terms. 

5) The rotor blade is modeled as a rigid blade with orthogonal springs 
located at the root of the blade (Fig. 3). This model enables one to 
represent simultaneously configurations employing either hingeless 

or articulated rotor system. The hinge sequence is given in Ref. 1. 

6) Since the geometrical nonlinearities due to moderate deflections of the 
blade are known to have significant role in rotary wing aeroelasticity 

The flexible portion consists of the elements having a length and shown 

in Fig. 2 . 
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[Ref. 2], these nonlinearities are included in the analysis. Retention 

of the nonlinear terms is based upon an ordering scheme [Refs. 1 and 2]. 

The blade degrees of freedom, representing the blade slopes, are assigned 

an order 0(e), where 0.1 < e < 0.15. The rigid body degrees of freedom 

3/2 

of the vehicle are assumed to be of a slightly smaller magnitude 0(e ) 

and the elastic deformations of the supporting structure are of the order 

2 

magnitude 0(e ) . This assumption is quite important for obtaining 

equations which are manageable from an algebraic point of view. The order- 

2 

ing scheme consists of neglecting terms of the order 0(e ) when compared to 
unity, thus 1 + 0(e 2 ) - 1. 

The equations of motion for the model vehicle (Fig. 2) are nonlinear coupled 
differential equations and they represent the coupled rotors/vehicle dynamics 
in forward flight. These coupled equations are classified in three groups, each 
group representing an appropriate sub-system equations. They are: 

1) rotor blade equations of motion in flap, lead-lag and torsion, 

2) rigid body equations of motion of the complete vehicle, 

3) equations of motion of the flexible supporting structure. 

The main advantage, due to separating the equations into various groups, 

luuoiolo ui me: auxxxty lU ucax wxlu vax xuuo ui piuuxcmo, ouen a. o iouxqxcu 

rotor aeroelastic stability or coupled single rotor/fuselage stability, etc., 
in a convenient manner. The coupled equations have considerable versatility and 
they can be used to study a number of diverse problems which are listed 
below: 


1. Isolate^ rotor aeroelastic stability. 

2. Coupled single rotor/fuselage dynamics. 

3. Response to cyclic and collective pitch inputs. 

4. Response to higher harmonic control inputs. 
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5. Stability analysis of twin rotor system connected by a flexible 
structure . 

6. Dynamics of a Hybrid Heavy Lift Airship. 

Depending on the type of analysis desired, the equations are simplified and 
modified to obtain an appropriate solution. 

Because of the unique nature of the multirotor model (Fig. 2), the results 
of the stability analysis could not be compared with any other results available 
in the literature. But on the other hand, experimental results are available, 
in the literature, for the aeromechanical stability of a single rotor helicopter 
in ground resonance. Hence, solving this problem analytically will provide an 
opportunity to validate both the equations of motion for the coupled ro tor/vehicle 
system and also the method of solution. Therefore, two different types of pro- 
blems are solved using the analytical model, for the coupled rotor/vehicle dyn- 
amics, presented in Ref. 1. In the first case, the equations of motion are used 
to predict the aeromechanical stability of a single rotor helicopter in ground 
resonance. It was found that the analytical results are in good agreement with 
the experimental results indicating that both the equations of motion for coupled 
rotor/vehicle system and the method of solution are valid. 

In the second case, the stability of a model vehicle (Fig. 2) representing 
an HHLA in hover is analyzed. The total number of degrees of freedom for the 
model HHLA consisting of two four bladed rotors is 31 . Hence, the stability 
analysis yields a total of 62 eigenvalues corresponding to these 31 degrees of 
freedom. Based on a careful parametric study, various blade and vehicle modes 
have been identified. A physical interpretation of the eigenvalues is obtained 
from a systematic study of the eigenvalue variations as a consequence of the 
variations of the vehicle system parameters. Finally the coupling between various 
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blade modes and vehicle modes is identified. 


In this report, the method of solution, the relevant trim (equilibrium) 
and linearized stability equations for the two applications mentioned above 
are considered and explored in detail. 


2 . METHOD OF SOLUTION 


The equations of motion representing the dynamics of the coupled rotor/ 
vehicle system, presented in Ref. 1, can be used to obtain either the response 
or the stability of the vehicle. The method of solution depends on the 
type of problem being considered, i.e. whether a response or stability analysis 
is required. For a stability analysis, one must distinguish between the case of 
hover which is relatively simple and the case of forward flight which is much 
more complicated. In this section, the method of solution used for the aero- 
elastic stability of a multirotor vehicle in hover is presented. 

The equations of motion, for coupled rotor/vehicle problem, are usually 
nonlinear coupled differential equations with periodic coefficients. These 
differential equations can be either ordinary or partial depending on the type of 
model used for the representation of the blade. If the blade is modelled as a 
rigid blade with root springs, the resulting equations will be nonlinear ordinary 
differential equations. On the other hand, if the blade is modelled as a flexible 
beam, the final equations will be nonlinear partial differential equations. In 
this case, the partial differential equations are first transformed into ordinary 
differential equations using Galerkin ? s method. Thereafter, the method of solu- 
tion is the same, irrespective of the modelling of the blade. In the present case, 
because the blade is modelled as a rigid blade with root springs (Fig. 3), the 
equations of motion are nonlinear coupled ordinary differential equations with 
periodic coefficients. To obtain the stability of the vehicle the following 
procedure is used: 

1. Evaluation of the trim or equilibrium state. 

2. Linearization of the nonlinear ordinary differential equations about the 
equilibrium position (linearized equations will have periodic coefficients). 
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3 . 


Transformation of the linearized equations with periodic coefficients 
to linearized equations with constant coefficients, by applying multi- 
blade coordinate transformation. 

4. Evaluation of the eigenvalues of the linearized equations with constant 
coefficients to obtain the information on the stability of the system. 
These four steps can be separated into two stages of analysis, namely, ( i) a 
trim analysis intended to establish the nonlinear equilibrium position of the 
blade, and (ii) a stability analysis of the linearized perturbation equations 
about the equilbrium state. A description of these two analyses are pro- 
vided in the following sections. 

2.1 Trim or Equilibrium State Solution 

In the trim analysis, the force and moment equilibrium of the complete ve- 
hicle together with the moment equilibrium of the individual blade about its root 
in flap, lead-lag and torsion are satisf ied respectively. It is important to re- 
cognize that only the generalized coordinates representing the blade degrees of 
freedom will have a steady state value representing the equilibrium position. The 
generalized coordinates associated with the rigid body motions of the vehicle are 
essentially perturbational quantities and hence their equilibrium, or trim, values 
are identically zero. In deriving the equations of motion for the flexible sup- 

structure, it was assumed that the vibrations of the structure occur about 
a deflected equilibrium position. The determination of the equilibrium position 
of the supporting structure is unimportant in the case considered here, for the 
following reasons: (a) this equilibrium position is not going to affect the 

equilibrium values of the blade degrees of freedom, since the blade equations 
contain only the time derivatives of the degrees of freedom representing the 
elastic modes of the supporting structure. The physical reason for this 
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mathematical dependence is due to the fact that blade inertia and aerodynamic 
loads depend on the hub motion and not on the hub equilibrium position. The 
hub motion is related to the fuselage motion and the vibration of the supporting 
structure, and (b) the final linearized differential equations used for the 
stability analysis do not contain any term dependent on the static equilibrium 
of the supporting structure because only the perturbational blade inertia and 


aerodynamic loads excite the vehicle rigid body motion and the vibration of the 
supporting structure. Hence, the generalized coordinates for the vibration modes 
of the supporting structure are again perturbational quantities. 

t* Vi 

The k blade degrees of freedom can be written as 

Flap 

Lead Lag (2.1) 

(ft = 4> k0 + Acf> k (40 Torsion 


B k ' 6 k0 + 46 k (W 
Sc ' A k0 + AC k <W 


where B^q, 4> k0 are the steady state values and AB k> A^, A<J> k are the per- 

turba t ional quant it ies . 

Linearization of the equations is accomplished by substituting these expres- 
sions into the nonlinear coupled differential equations and neglecting terms con- 
taining the products or squares of the perturbational quantities. The remaining 
terms will have either the steady state quantities as coefficients or the time 
dependent perturbational quantities, multiplied by the steady state values 
or some appropriate constants. Separation of these terms yields two 
groups: one group of terms contains only the steady state quantities and con- 

stants (i.e., time independent quantities). These represent the trim or equilibrium 
equations. These are nonlinear algebraic equations which represent force and 
moment equilibrium equations determining the steady state. The second group con- 
tains the time dependent perturbational quantities and represents the linearized 
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equations of motion about the equilibrium position. These linearized dynamic 
equations of equilibrium are used for the stability analysis. The steady state 
equilibrium equations can be written symbolically as 
for the complete vehicle 

F = F = F = 0 (2.2) 

x v z 

M =M =M =0 (2.3) 

x y z 

and for the individual blade 


M =M =M , = 0 

B C <f> 


(2.4) 


where F , F and F represent the forces of the vehicle in X.Y.Z directions, 
x’ y z 99 9 

respectively; represent the moments on the vehicle about X,Y,Z axes res- 
pectively; and represent the moments of the blade forces about the root; 

respectively . 

In these equations, F^, F^ and are identically zero. The remaining 


equations for the vehicle can be written as 


F = T, + T 0 + P3 - W = 0 


z 

M = 0 

y 

M = 0 


(2.5) 

( 2 . 6 ) 

(2.7) 


where T^ and T^ are the magnitude of the thrust developed by the two rotor systems 
S 

R^ and is t ^ ie static buoyancy on the envelope, and W is the weight of 

the complete vehicle. 

The quantities T^ and T^ are functions of the steady state flap, lag and torsion 
angles, collective pitch angles and the operating conditions of the rotors. 
Equation (2.7) for M z represents the torques developed by the two rotor systems. 
These torques can be either balanced by having a tail rotor for each main rotor 
or by having two counter-rotating main rotors. For the case where the rotors 
are assumed to be counter-rotating, the blade loads are to be evaluated for 


the two rotor system with angular velocities + fie^ and - fie^ respectively. 
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In Ref. 1, the rotor loads are derived for a typical rotor with angular velocity 
+ fte z and the same expressions for loads are used both rotors. Thus we assume 
that the torque developed by each main rotor is compensated by a tail rotor. 

Equation (2.6) for consists of the pitching moments developed by the 
thrust of the rotors and gravity loads on the various components. 

The steady state moment equilibrium equations for the individual blade 
will have the following symbolic form 


M e 

-4 


r 1 

hcO’ 


Q± o> 

= 0 

(2.8) 

M c 

■4 


r* 

ScO’ 

'ho’ 

e o ) 

= 0 

(2.9) 

M 4> 

-4 


ScO’ 


e o> 

= 0 

(2.10) 


where i = 1,2 refers to the two rotor systems and R 2 respectively 
and k refers to the k ^ blade in the i th system. For the case of steady 
state, all the blades in each rotor system will have the same steady state 

values (or equilibrium quantities) ana thus tne suoscript: ! k ! can be deleted. 

Equations (2.5), (2.6), (2.8) - (2.10) are nonlinear algebraic equations. 
There is a total of 8 equations and 8 variables (Sq> Cq, 0q ; i = 1,2). 

These eight equations can be solved iteratively by the Newt on-Raph son method, 

to obtain the steady state values. Failure to converge during iteration can be 
attributed to the divergence or static instability of the blade. (These equili- 
brium equations are given in the next chapter.) 

In deriving the equations of motion, the inflow ratio X is assumed to be 
constant over the disc. The typical value chosen for the inflow ratio is its 
value at 75% of the blade span. It is given as [Ref. 3] 



where 


0q is the collective pitch of the blade. 


( 2 . 11 ) 
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2.2 Iterative Procedure for the Trim Solution 


The equilibrium equations of the blade (Eq. 2.8 - 2.10) and the equations 
of the complete vehicle (Eq. 2.5 and 2.6) have to be solved numerically to obtain 
the steady state values of the blade deflections in flap, lag and torsion 
(Bq, Cq, 4>q ; i " 1,2) and the collective pitch angles ( 0 q, i = 1,2) of the 
rotor systems and . Blade equilibrium is obtained by an iterative procedure. 
It can be seen from the equilibrium equations of the blade (Eq. 2.8-2.10) that the 
equations for the blade in i ^ rotor system consists only the variables corres- 
ponding to that rotor system. Hence these equations can be solved separately for 
each rotor. By assuming a collective pitch angle of 0 q, the equations 

(2.8) - (2.10) are solved to obtain the equilibrium angles (Bq, Cq, 4>q) of the 

2 

blade in rotor system R^ . Then by assuming a collective pitch angle 0^, the 

2 2 2 

equilibrium equations are solved again for the equilibrium angles (Bq, Cq, (^q) 
of the blade in rotor system R^ * It is important to recognize that these equili- 
brium angles of the blades can be different for the two rotor systems R^ and R^, 
because these angles depend upon the operating conditions and the blade parameters 
which can be different for the two rotor systems. However all the blades in each 

rotor system will have the same equilibrium angles. After obtaining Bq, £q, $q 

- 22 . 2 . 

and Bq, Cq, these quantities are substituted in the vehicle equilibrium 

equations (Eq. 2.5 and 2.6) and these two equations are solved simultaneously 

1 2 

to get the updated values for the collective pitch angles (0 q and 0 q) for the two 

rotor systems. With these updated collective pitch angles, the blade equilibrium 

equations, for the two rotor systems, are solved again to obtain a new set of 

equilibrium angles for the blade (Bq, Cq, ^q ; 1 = 1,2). These equilibrium angles 

of the blades are again substituted in the vehicle equilibrium equations to get 

1 2 

the second stage updated values for the collective pitch angles 0 q and 0 q. These 
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steps are repeated until convergence is achieved. This procedure is also il- 
lustrated by the following flow chart. A computer program implementing this 
calculation for the trim (or equilibrium) position of the blade was developed. 

A check on the number of iterations. is provided in the computer program to avoid 
excessive use of computer time in case of divergence. 


- i = 1 

7 . 1 

■ Assume 0^ - 

i. u » 

t 

Evaluate A . 

i 

^ using Eq . (2.11) 


Solve equations (2.8)— (2.10) 
iteratively for 

r 


If i < 2 
i = i + 1 


True 


False 


Using Bq, Cq, 4>q J 1 = U 2 
solve equations (2. 5) -(2. 6) 
for 0 q and 0^ 


False 


Check for convergence. Convergence is 

achieved if the differences between 

12 -5 

initial and final values of 0 q and 0^ are < 1 .0 x 10 

k--' - t • *~ 

True 


Stop 


Flow Chart Illustrating the Solution Procedure for 
the Blade Equilibrium Position 
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The procedure described above is restricted to the case of hover where the 
coordinates describing blade equilibrium are not time dependent. For the case of 
forward flight, the equilibrium values will be time dependent and a more complicated 
procedure, described in Ref. 4, is required to determine the trim quantities. The 
basic difference is that, for the case of hover the trim values can be obtained by 
solving a system of nonlinear algebraic equations, while for forward flight the 
solution of a system of nonlinear coupled ordinary differential equations is 
requ ired . 

2.3 Description of the Stability Analysis 

The perturbat ional equations of motion, linearized about the equilibrium 
position, can be written in the following form 

[M] {q} + [cl {q> + [K] {q> = 0 (2.12) 

where {q} contains all the degrees of freedom representing the blade motion, 
the rigid body motions of the vehicle and the flexible modes of the supporting 
structure. 

The matrices [M] , [ C] . [K] can be identified as mass, damping and stiffness 

matrices respectively and the elements of these matrices are functions of the 
equilibrium values. 

The stability of the vehicle about the trim condition is obtained by solving 
the eigenvalue problem represented by Equation (2.12). For convenience Equation 
(2.12) is written in state variable form 

{y} = [F] {y} (2.13) 

where (y} T = 

and (y^ } = {q} ; { y = {q} 
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and 

-IM]" 1 [C] 

i -i 1 

| - [M] X [K] 



[F] - 

|— 



[I] 

1 0 
I 


Assuming 

a solution of the form {y} = {y} 

c i> 

e , Equation (2.13) reduces to a 

standard 

eigenvalue problem 




[F] {y} = S{y } 


(2.14) 


The eigenvalues of Eq. (2.14) can be either real or complex conjugate pairs. 

s k-°k liw k <2 - 15> 

The complex part of the eigenvalue (a)^) represents the modal frequency and the 
real part (cr^) represents the modal damping. The system is stable when < 0 
and the stability boundary is given by = 0. 

This relatively simple procedure can become complicated depending on the 
form of the matrices [M] , [C] and [K] • In the aeroelastic stability analysis 
of isolated rotor in hover, these matrices contain constant elements. Thus 
solution of this eigenvalue problem is straight-forward. However, in the case of 
coupled rotor/body system stability analysis in hover or for stability of iso- 
lated rotors in forward flight, these matrices will have elements which are time 
dependent. The reason for the appearance of time dependent or periodic coef- 
ficients, for these two cases, is different. For coupled rotor/body problem, 
these matrices become time dependent due to the fuselage perturbational motion. 
This fuselage perturbational motion introduces, through the hub motion, periodic 
terms in inertia and aerodynamic loads of the blade. In the case of isolated 
rotor in forward flight, these matrices become time dependent due to the per- 
iodic aerodynamic excitation associated with forward flight. 

When the coefficient matrices of the linearized perturbational equations 
are periodic the stability analysis can be performed by applying 
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one of two possible techniques. One can use either Floquet theory or introduce 
a multiblade coordinate transformation [Refs. 3,5]. For the coupled rotor/ 
body type of analysis in hover, which is the main objective of this study, 
the multiblade coordinate transformation is successful in eliminating the 
time dependent coefficients from the equations of motion. During this coordin- 
ate transformation, the blade degrees of freedom in the rotating coordinate 
system are transformed into a nonrotating hub fixed coordinate system. It is 
worthwhile mentioning that this transformation is also frequently denoted by 
the term Fourier coordinate transf ormat ion , Coleman transformation and more 
recently rotor plane coordinate transformation. 

The multiblade coordinate transformation is implemented by applying the oper- 
ators, given below, to the blade equations. 


1 

N 

N 



E (•••) 


collective operator 

k=l 




N V 



1 

£ <-D 

.) 

alternating operator 

N 

k= 1 



1 

N 

N 



E cos nip 

(...) 

n-cosine operator 

k=l K 



1 

N 

N 



E sin nip, 
k=l 

(...) 

n-sine operator 

where N is the number 

of blades 



and n = 

1 , . . . . L 

L = 

2 for odd N 



L = 

N " 2 * 

2 for even N 


The resulting equations are identified according to the operator used in the 
transformation. These operators are applied only to the blade equations because 
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only the blade equations are written in the rotating coordinate system. The 
vehicle equations are written in the nonrotating frame and in these equations, 
rotor loads appear as a sum of various loads due to the individual blade. A 
clear description of the theory of this transformation is presented in Refs. 3 
and 5. A brief summary of this transformation is also given in Appendix A. 

During the derivation of the final trim and linearized perturbation equations, 
two items are noteworthy. First, in Ref. 1, the blade loads are 
derived for typical rotor blade rotating with angular speed ft. During 
nondimensionalization of various quantities, the time is nondimensional ized as 
= ftt . When the general expressions for the blade loads are applied to two 
different rotors operating at different values of ft, then the nondimensional time 
will be different for the two rotor systems. Consider the two rotor systems 
R^ and R 2 to be operating at angular speeds ft^ and ft^ respectively. Then ft^t 
is the nondimensional time used in the rotor load expressions for the rotor system 
R^ and ft^t is the corresponding nondimensional time for rotor system . For 
the sake of consistency, the nondimensional time should be made the same for all 
rotor systems. Assuming that ft^ is the reference R.P.M. for the non- 
dimensionalization of time. Then the time derivative termswhich appear in the 
blade loads of rotor system R^ must be multiplied by a factor ( — ) . The power 
of this factor depends on the order of the time derivative. 

A similar problem is also encountered when multiblade coordinate 
transformation is applied to a multirotor system. In the n-cosine and n-sine 
transformation Eq. (2.16), ij^ refers to the azimuth angle of the k t ^ 1 blade 

i) k = ftt 4- 27Tk/N (2.17) 

where it is understood that ft is the angular speed or R.P.M. corresponding to 
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the particular rotor, and N is the number of blades in the rotor system. 

In the case of multiple rotor systems, this equation can be written for the 
k t ^ 1 blade in the i ^ rotor system as 


. i n , 2uk 

\ * V + 


(2.18) 


th 

Note that the above equation contains the R.P.M. or angular speed of the i 

rotor system. This value of Eq. (2.18), must be used in the n-cosine and 

th 

n-sine transformation operators for the i rotor system, and care should be 
exercised when applying these operators in transforming the time derivatives 
of the blade degrees of freedom. For consistency in nondimens ional izat ion , 
the time derivative terms in the i^ rotor system are nondimensionalized with 
respect to the reference angular speed of the rotor . These statements 
imply that for multirotor systems the expressions provided in Ref. 3 and 5, 
for transforming the time derivatives of the rotating blade degrees of freedom 
to the nonrotating system, should not be used directly. The correct form for 
implementing this transformation for a multirotor vehicle where each rotor is 
operating at different angular speeds, i.e. and provided in Appendix B. 

Consistency in nondimens ional izat ion of time and the mulliblade coordinate trans- 
formation for multiple rotor systems can be both achieved by multiplying the first 
and second time derivative terms, in the transformed multiblade coordinates for 


n 


2 


.th 


1 1 

the i rotor system by ( ) and ( ) respectively, where Q ± is the i 

i i 

rotor R.P.M. and is the reference angular speed. 

(2) The second noteworthy item is related to the rotor hub loads. When deriv- 
ing the equations of motion of the vehicle, the rotor hub loads have to be eval- 
uated. The rotor hub loads are obtained by summing up the contributions from the 
individual blade loads. The expression for the individual blade load will have 
the centrifugal term as the leading term, the order of magnitude of this 
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term is 0(1) . After summing up the individual blade loads, the resulting hub 

load expressions will have a leading term of order of magnitude 0(e) only, 

because the centrifugal contributions from the blades cancel each other and 

thus the net contribution due to these terms is zero. Therefore, care must 

5/2 

be taken to retain terms up to order 0 ( £ ) in the individual blade load ex- 

pressions, so that the resulting coupled rotor/vehicle equations of motion will 
represent a consistent nonlinear mathematical model. 

The two items, discussed above, have been carefully implemented in the 
equations which have been derived in this report. Next these equations are 
specialized to study air resonance type problems. For this class of problems, 
it is common practice to suppress the vertical motion motion and the yaw 
degree of freedom, thus in the final equations 


and 


R - 0 , R = 0 , R =0 
zs zs zs 


e = 0,0 = o ,0 =0 


zs 


zs 


z s 


(2.19) 

( 2 . 20 ) 


have been substituted. The final equations for the equilibrium position (trim) 
and for the stability analysis are given in the next chapter. 
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3. EQUATIONS FOR TRIM AND STABILITY ANALYSIS 


In this chapter, the complete set of equations, used for the equilibrium 
position and stability analysis, are presented. The linearized stability 
equations are given in multiblade or rotor-plane coordinate system. This 
chapter is divided into two major sections. In the first section, the equa- 
tions pertaining to the twin rotor model (Fig. 2) representing an HHLA, are 
presented. Two sets of equations are provided for the HHLA model: one for 

articulated rotors and the other for hingeless rotors. The second section 
presents the equations used for predicting the aeromechanical stability of 
a single rotor helicopter in ground resonance, including the effect of the 
aerodynamic loads. 

3.1 Equations for Twin Rotor Model of an HHLA 

The degrees of freedom included in the analysis of twin rotor model of 
an HHLA are listed below. 

Blade degrees of freedom 


Flap 


4 * 

4c • 

3 1 

ns 


Lead-lag 

4 • 

^-M ’ 

4* 

s 1 

ns 

i = 1,2 refers 

Torsion 

4) 1 , 

*4 • 

4c ’ 

4) 1 

ns 

systems 


The subscript M refers to the collective mode, -M refers to the alternating 
mode (only for rotors with even number of blades) , nc refers to the n-cosine 


mode, ns refers to the n-sine mode. 

Rigid body degrees of freedom of the vehicle 


X - Translation 

R 


xs 

Y - Translation 

R 


ys 

Roll 

9 


X 

Pitch 

0 


y 
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Degrees of freedom for the flexible supporting structure 

r 

Bending in X-Y plane (Horizontal plane) 

Bending in X-Z plane (Vertical plane) 

Torsion F 

o 

These degrees of freedom represent the three normal modes of vibration of the 
supporting structure, two for bending in two orthogonal planes and one for torsion. 

Thus the total number of degrees of freedom for a vehicle consisting of 
two four bladed rotor systems is 31. 

3.1.1 Static Equilibrium Equations (Trim Equations) 

In the following, the nonlinear algebraic equations required for the calculatior 
of the trim quantities, for a hovering vehicle, are presented. These equations are 
solved to obtain , 6 q and A where i=l,2, refers to the 

two rotor systems. Since the form of the equations are the same for the blades 
in both rotor systems and , only one set of blade equations is needed. 

It should be noted that these blade equations are solved separately for each rotor 
system with its own parameters. Furthermore for convenience in writing these equa- 
tions the superscript , i, is deleted from the equilibrium quantities only. 

Flap Equation 

^kO + ^ko + ^ko F 'd 3,i - > 

+ 6 k0 hcO F T (4,l) + ^kO ^kO F T (5,l) + ^kO 6 k0 F T + F x = 0 

(3.1) 

i = 1,2, refers to the two rotor systems and 

-3 -2 

-2 -2 -2 2 0 Q 

F T (l,i) = 0) F + (co L - (Op) sin 0 Q + -y + e y 

—2 —2 0*^ 

F (2 , i) = (a) - co ) sin0 cos0 + V —r~ £3 

i L F 0 0 4 p 
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F T (3,i) = 

£ 4 l 3 

" v ( X + T 

2e ) 

F x ( 4 ,i) = 

£ 4 

V T 


F x (5,i) = 

-2 -2 -2 

- W F " ( \ “ V 

sin °o 

F T (6,i) = 

-2 -2 

(o> L - (Op) sin0Q 

COS0Q 

F x (7,i) = 

£ 3 £ 2 

V V + !T e > 

r ^ 

- v [ X 


where 


-2 

W F 


K, 


2 3 
mJTR 


-2 

W L 



2 3 
R 


£ = 1 - e 


£ 3 

9 0 + T 


(-X+2 e 0 q) 



iX ] 


e 

V 

Lead-Lag 


_e 

R 

p . abR 
A 

m 

Equation 


4o 4(1,1) + 6 k0 4(2,1) + ^ k0 4(3,1) 


+ ^kO 6 k 0 

L x (4,i) 

+ <P 

O 

M 

k_p 

O 

L x (5,i) 

+ 6 k 0 ? k 0 

L x ( 6 ,i) + 

V 7 

,i) = 0 









(3.2) 

L x (l,i) = 

- “£ + 

(5 L 

-&■ 

sin20 o ' 

l 2 - 

~2 e + v(- 

£ 4 

■ 7 - $ e n - 
4 p 0 

l 3 

3 

2 x e ) 

p 

L x ( 2 ,i) = 


$ 

sin 0 ^ 

cos0q 





L T (3,i) = 

, £ 3 
v (- T 

X - 

2 X 

e) 
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V 


L T (4,i) 


-2 -2 -2 2 
- 0J L + (co - co F > sin 0 q 


L T (5,i) 


-2 -2 

(y - a> ) sin0Q cosQq - 


t 


L T (6,i) 


£ 3 i 

v y x 


L T (7,i) 


c A — 3 

v{ - ( -7~ + 2 -r- e ) 

a 4 3 


f>V T 


Torsion Equation 


4» k0 T x (l,i) + £ vn T t (2 , i) + 3 vn T t ( 3, i) 


= kO T 


kO T 


kO 

C k0 

T x (4,i) + 

^kO ^ 

kO 

r T ( 5 

,i) 

+ ^kO 

T x (6,i) 

kO 

4 

T x (7,i) + 

4 5 

3 , 
kO 

r x (8 

> 1) 

+ \o 

‘’kO T t ( 9,i) 

2 

kO 

? k0 

V 10 ' 1) + 4 

^kO 

^kO 

T 

T 

(11, i) 

+ 6 k0 C k0 ^kO 

'kO 

^kO 

\o V 13 ' 

( i) + 

B k0 


\o T T (14 - i) + 6 k0 *k0 

’kO 

4o 

♦m V 16 ' 

>i) + 


(17, 

i) 

= 0 



,, -2 -2 - , . 2 A 2„ w X MB3 I MB2 

T (1, i) = - 0) - <o + l (sin e Q - cos e Q ) ( — y - — y 

l l l 2 mR mR 


-2 “3 -2 

+ x sin0 hr B + Vx ( h- + hr 2 e ) 
I 0 2 p A3 2 


p3 n2 n 2 a-4 _ _ 

T T (2,i) = “ B ( V + V 5 ) - *1 Sin9 0 V + V( V 6 0 + V ( ^ +2ee 0 } " T eX) 


* 3 


(3.3) 


a 2 


- , £ 3 „ x , 77 , X MB3 2 A . X MB2 . 2 Q x 0 

+ Vx A (- y 6 p ) + A ( yy cos 0 q + — y sm 0 Q ) B p 


mR 


mR 


c -4 -3 

C d0 , Jl . , l' 


i 3 


i 2 


_ £ 2 


T T (3,i) = - V[- -f - ( y + 2 5) - y A0 O - i- X(-X+e0 Q ) ] - y y 


cosG, 


m ,, -2 9? . T/ I MB3 2 a . X MB2 . 2 n . , ,, l* 

T t (4,i) = 0> T - y + £( — — cos 0 Q + — y sm 0 Q ) + V y 

1 mR mR 


P 0 
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t ! - I 3 

- V T 2As p - vx aT 

-2 -4 -3 

T T (5,i) = - x T cosQq 4y + v ( y- + y- 2 e ) 

T T (6,« = -v6 p ^ 

I 4 


c d0 7 £ 4 . „ £ 3 - 


£ 3 


£ 2 


T x (7,i) = -V— - v[ - — ( y- + 2 i ) - A e Q - A(-A+e0 n ) ] 

r n . n r 


7t4 


■?r3 


T x (8,i) = - v(- y- 0 Q + 2 ^ A ) 


t 


£ 3 


T x (9,i) = V- B p e 0 - vy 2 A 8 p 

t l 3 

T T (10,i) - - V (- y- 0 Q + 2 V A) 


T x (ll,i) = V^- 

T x (12,i) = v^- 
£ 4 

T x (13,i) = vy- B p 

V 14 *« - 6 p 

-3 -2 

T x (15,i) = v(yHyAi) 

-3 -2 

T x (16,i) = v(y Uy A i ) 

-2 II 

T x (17,i) = X x cosBq [- B ( y- + l S) ] + £ [-( -yp- - -yy ) sin0 o cosQy] 

^ mR mR 

-3 -2 

+ ( 3 0 O + A+ 2 e 0 Q ) - i £A ) 

Equations (3 . 1) -(3.3)are valid for both articulated as well as hingeless blades. 

For articulated rotors ol = 0 

F 


oo =0 

L 
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0 


and 


For hingeless rotors 


and 




n 2 3 
K, 

4> 

2 3 
mfl R 

0 



K . Ka 

♦b 

K, +k a 

4> b 4> 


c 


In the above blade equilibrium equations, the inflow ratio X is [Ref. 3] 


„ / 1.50 n 

= it < - 1+ V 1 + wify > 


(3.4) 


Force Equilibrium Relation 

For hover, the thrust developed by the two rotors and the buoyancy force 
on the envelope must balance the weight of the complete vehicle. Also, the 
pitching moment due to the various forces about Y-axis (Fig, 2) must be 
zero. 


£ [ £ mf2R {v[ 3 (0 Q + ^ k o ^kO B p " 5 k0 6 k0 ) 

1=1 k=l 

l Z - - - - , 

+ i-(-A+2e0 o +2e <t> k() ) - l e X ] } ] . 


+ P Z ( W F1 + W F2 + W UN + W EN + ° 


(3.5) 


Moment Equilibrium Relation 


£ I- < £ « 2r3 <v l -V <e„ + - c k0 \o> 


i=l F± k=l 


a l 2 


kO ^kO M p ^kO M k0 


(-X+2e0 o + 2e <f> kQ ) - l e X ] > } ] 


+ *F1 W F1 + *F2 W F2 + h 5 W s = ° 


(3.6) 
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In equations (3.5) and (3.6) the symbol 1 i f outside the brackets indicates 
that all the quantities within the bracket refer to the i fc ^ rotor. The 
quantities within the bracket can be different for different rotor systems. 

In total, these are eight equations, three blade equilibrium equations 
for each rotor (total 6) and two vehicle equilibrium equations. As explained 
in the previous chapter , these equations are solved iteratively to obtain 

the equilibrium state. The solution consists of , 

0^ ; 1-1,2, these equilbrium quantities are the same for all the blades 

in one rotor system. 

3.1.2 Stability Equations 

The equations of motion for the blade as well as for the vehicle are 
linearized about the equilibrium state. These linearized equations are then 
transformed into multiblade or rotor plane coordinates. The final linearized 
equations, written in the multiblade coordinates are given below. 

Collective Flap Equation 


4 *.<!.« + 4 F c <2 ’ 1) + 4 F c (3 ' i) 


+ 4 + 4 F c (5 - i) + 4 F c (6 ’ i) 

c 

+ B* F c (7,i) + e y F c (8 , i) + F c (9,i) 


+ 0 F (10, i) + F (11, i) = 0 

y c ’ R c v * 


(3.7) 


th 

i = 1,2 refers to the i rotor 


F c (l,i) = Up + (uj^ - u^) sin 2 e Q + <t> kQ (a^ - w 2 ) sin6 Q cose Q 

r . . li - 

+ V 4 5 k0 + 3 + 2 e 


F c (2,i) 


—2 —2 —2 2 —2 —2 
- 4> k0 ( Wp + ((-^ - Wp) sin 9 q ) + (a^ - aip) sin0 Q cos0 q 
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+ v x <8 p + f W 

F c ( 3,l) = - ; k0 { ^ + (^ - ^> sin 2 e 0 ) + B k0 (Sj - “f> si " e o cos0 O 

I 3 , - <i 4 

- V -=— 2 e - V ~r~ 

3 4 

£ 4 I 3 n l 

F c (4 - l) - (S SF + v x + V— e) jj- 


$3 £ 3 

F c (5,i) = [2 — (6 k0 + B p ) - 2v X (9 0 + \o ) + ^7 

£ 3 _ i _ 

F c (6,i) =[- vyb-yvb — cosQ Q ] jj- 

o 3 i J 3 _ 2 

F e (7>i) = [ 3 + I V T b COS0 O 1 ( ST } 

F c ( 8 ,i) - - * pi V- ( n- ) 

1 

0*3 SI . 1 

F c (9,i) - n 2 « rl ) v T x 


F c a°,i) - - Vi T < 57 > 

1 

,2 Q, 2 

f (11,1) - n 2 (» F1 ) - ( 57 ) 

1 

It should be mentioned that the angular speed SI can be different for the two 
rotors, and thus the blade static equilibrium represented by 3 kQ , C kQ , <t> kQ , 

0q, A can differ from one rotor to another. Thus the coefficients can 
be different for the two rotors. Actually Eq. (3.7) represents two equations, 
one for each rotor system. 

Alternating Flap Equation (For even N only) 

®-M F A (1>1) + ’-M F A (2 ’ W + *-M V 3 -” 

+ F A (4,i) + I (5,1) + F A (6,i) 


+ i.‘ F A (7,i) - 0 


(3.8) 
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F A^ » ^ = + (<\ ~ Wp) sin 2 e 0 + 4> k0 (to 2 - w 2 ) sin0Q cos0 Q 


£ 4 £ 3 l 2 - 

+ v_ 4 ~ ? kO + T + T e 

_2 _2 -2 2 —2 —2 
V 2 ’i> = “ ^kO ^ W F + “ W F^ Sin 6 o } + K ” W F^ sin0 o cos0 O 

+ v x « P + V 

F A ( 3,i) - - c k0 < “p + <\ - “p> sin 2 e 0 }+ e k0 <(7ij - ^)sin6 o cos0 o 

- vA-2e-v^- 

- if* 7 2 

F A (4 ’ i} = (8 sf +V X + V T e) ^ 

F A (5,i) - { 2 £ ($ k0 *( f )-Zv(( V * t0 ) + V £■ X }^i 

o 3 _ 1 _ ?3 Q X 

F A (6,i) = { - v— b - y V b T COS0 O } JT 

i 

V 7 ’^ = { T + T cos6 o } ( TT> 

i 

n-Cosine Flap Equation 

3* F (1, i) + 6 1 F (2, i) + C 1 F (3, i) + ^ F (4,i) 
nc nc ns nc ’ ^nc nc 5 b ns nc 9 


t \ J_ C TT ( Q 


ns nc 5 


i i 

-U (f) 17 i'i X T? _L TT ( 7 

_ *■ v . ’T A \D * - l / 1 KJ -l 

nc nc ns nc nc nc 

+ *4 F „c (9 > 1) + 4 F -< 10 - I > + 4 F„„(u,i) + e„ F„„<12,i) 


nc nc 


nc nc 


: c, 

+ e „ F ^< 13 -« + e x F „ c < 14 .« + R „. F „( 15 .« + -r F „<i6,i) 


y nc 


ys nc 


y nc 

fi 

R nc 


+ C 3 F nc (l 7> i) = 0 


(3.9) 


F (1, i) 
nc 


-2 . ,-2 -2 


-2 -2 


= U)p + (w T - W p ) sin 0 n + ^(dl, - ) sin0 n cos0 


0 Y kO L 


'0 — 0 


-4 -3 -2 -3 -3 

. „ i _ . a - 3 . a - 2 sl 21 ,-r 

+ v ~ c ko + T + T e _n T “ n 2 vb l" cos9 0 
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t l 2 - - s 

F_(2,i) = n( V - 7 - + V — e + g sp ) 


nc 


F nc (3,i) = - <t> k0 { “f + ( ^L ' “f 5 Sln ^ 0 O } + (W L " W F } Sin0 O COS9 0 


-2 - 2 , 


tt3 


+ v (B p + e k0 ) 
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1? 


F nc (4,D = n { 2 f- (B k0 + B p ) - 2 v (6 0 + ^ + v X } 


,-2 - 2 , 


F nc (5,i) = - C k0 ( ^ + (w 2 - sin 2 0 o } + B k0 (o£ - Up sine Q cos0 o 
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? 3 _ l l 3 

F nc (6,i) = n {- V -y- b - j V b — cosSq } 
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? 3 i Z 3 - 
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F nc ai.i) = { T + \ V T E cos0 o } ( sT } 


l 3 2 

F nc (12 ’ 1) = “ 6 n T ( ft7> 


£ 4 fi l 


F (13, i) = - 6 V -j- q- 
nc ’ n 4 ft. 


-3 -3 ^2 ft^ 

F nc <W.« - < + «„ 2 T + S n 2v T 9 0 E 2- S 2 V - X6 n ) ^ 

i 3 J 2 ^1 

F nc (15,i) = (- 6 n 2V^- e Q + VyH n ) ^ 
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F nc (16 ’ i) = ( - 6 n 2 V T 6 0 W + *1 a Fi } V T X V C 
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*3 


*2 


n. 


j? - 1 

F nc (17,l) = ( 6 n 2V T 9 0 h 2 n 3 (J W " h 2 n 3 ( *Fi* V T“ X V IT 


where 


6 = 1 

n 

- 0 


when n = 1 
n ^ 1 


n-Sine Flap Equation 

B* F ( 1 , i) + B* F(2,i) + ^ F (3,i) + C 1 F (4, i) 
ns ns nc ns ns ns nc ns 

• • • ♦ # 0 

+ <J> F (5, i) + (j) 1 F (6, i) + B 1 F ( 7 , i) + B 1 F (8,i) 
ns ns nc ns ns ns > nc ns * 

• *• 

+ iL F „.«.« + F ns < l °.« + eh F „.ai,i) + e. F„(i2,i) 


ns ns 


x ns 


U 


+ 9 x ? n S W ’» + 9 y F n S <l4 -« + \s F n S (15 ' 1) + f F ns (16 - 1) * 0 < 3 - 10 > 

F ns a,i) - “f + ( \ " si “ 20 o t +k0 ( \ " “f 5 sin9 0 Cos9 0 

-4 -3 -9 -3 -3 

. „a _ . a . a - 2 r i c r 2 

+ v X ? ko + T + T e " n T -i vb T n cose " 

A 4 I 3 - 

F ns ( 2 » i} = n { - - vf- e - g gF } 

F ns (3,i) = “ ^kO { + (i \ ” ^F } sin20 O } + ( “L _ “f } Sin9 0 C ° S0 O 

A 4 


£ 3 


+ v 1- (6 P + W 
£ 4 


£ 3 


F ns ( 4 ,i) = «{ - 2 T (B k0 + 6 p ) + 2 V (6 q + 4> k0 ) -V-j- A } 


F ns^ 5>1 ^ “ ^kO ^ “f + ^ Sin 6 0 } + ^kO ^ “ W F^ sin0 O cos0 O 

A 4 £ 3 9 - 

-v-7- - v - 5 - 2 e 

4 3 

-3 -3 

F ( 6 , i) = n { v\- b+|vb 4 - cos6 n } 

ns 3 2 3 0 

jf* n3 _ 

F ns^ 7 , "4* + V T 6 + g SF * TT 

i 


29 


9? 1 _ 

W-W ■ n( - 2 T -l 2vb T cos9 0 > 

F ns (9 -« ■ ' 2 T (6 kO + V - 2v T < 0 O + +kO> + v T x >ST 

1 


4 v "0 ’ Y kO / 


0 ^ _ ] _ ? 3 ^1 

F ns (10,i) = { ~ V T b - 2 v b T COS0 O } IT. 

X 


J 2 1 - £ 2 2 

F „s< u - 1 > ■ < T + 1 v b T “"V ( st > 

1 

1? 2 
i 

^4 

F n S (13 - i > ■ s » v T<aT) 

F „s <14 ’ i > - { + 6 „ 2 T + «n 2V T e 0 K 2 - K 2 v T X6 „ > ^ 

£ 3 I 2 

F (15,i) = (6 n 2v^ e.-VyH) ^ 
ns n 3 0 2 n 

_3 _2 

F ns (16,l) = R{ - <S n 2 v— 0 O h 2 n 2>x (^ F± ) + h 2 n 2>x (^ ± ) V— X6 n > 
Collective Lead-Lag Equation 


4 L c (i ’ i} + 4 L c (2 » i} + 4 L c (3>i) 

+ 4 v 4 ’*) + 4 L c (5 » i) + 4 L c (6 ’ i} 

. ... • E, 

+ 4 L c (7 » i) + 4 L c (8>i) + 6 y L c (9ji) + ~R V 10 * 15 = 0 O »11) 

L c (l,i) = - w 2 + (4 - 4) sin 2 0 o + <f> kQ (4 - 4> sin0Q cosGq 

l 2 - £ 4 

. _ e . v _ e o g p 

L c (2,i) = - <J> k0 { 4 ~ (4 - 4> sin 2 6 0 } - sin6 0 cos9 0 


L c (3,i) " ■ e ko { 4 ■ ( 4 ■ 4 } sin 20 o }+ \o ( 4 ■ 4 } sin6 o cos9 o 

* 3 , 

- V-y \ 
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C dO 0 ^ - ^1 

L c (4,i) - (- 2v -^ A - - vi- e 0 x - g SL ) jj - 


L c <5,:t> - {2 4 <B k0 + 6 ) - v-£ <e 0 + * k0 ) - (-2X + 5 6 0 ) } Ai 


1 - 2 3 
L (6,i) = - v b — sin0 


o n. 


j 3 n 2 

V 7 - 13 ■ - T < al > 

1 - I 3 n l 2 

L (8,i) - 2 vb T Sin0 o ( JJ- ) 

1 




fi. 


L„(9,i) = - Kt { - V 4- (e 0 +4> k0 ) + v hr 2Ur 


Fi 


I 3 


£ 2 


fi . 


L c (10,i) = n 2 (^ Fi ) { -V^J- (0 O + 4> k0 ) + V \ 2 X } 

Alternating Lead -Lag Equation (for even N only) 


w . 


L A (l,i) + 6_jj L A (2,i) + 4>_i L a (3 , i) 


’-M A' 


-M “A' 


?-M L A (4 - 1) + L A <5 -« + +-M "A* 6 - 13 


+ c.‘ V 7 - 13 + b -m V 8 ’ 13 ■ ° 


(3.12) 


L A (l,i) = - U) 2 + ( U) 2 - W 2 ) sin 2 0 o + 4> k0 (Si£ - Up) sin0 Q COS0Q 

£ 2 - 7^ 

- T e - v T Vo 


L A ( 2 , i ) 


" \o * “l “ ■ “jp sin20 O * " “ W F^ Sin9 0 cos9 0 


L A (3,i) 


-2 ,-2 - 2 » . 2 


-2 - 2 , 


* B k0 { " (W L " W F ) Sin 6 0 } + ^kO ( ^L " V Sin9 0 COS0 


l 3 

- V- A 


L A (4,i) 


L A (5,i) 


= { - 2 v 


£ 3 


c d 0 £ 4 


£ 3 


- 


Xj ^ r \ y “ i a 

/. — ^ o ^ Sot ■* O 


3 0 

£ 4 


SL fi. 
x 

£ 3 


fi 


1 


^ 2 3 ^ 0 kO + V V 4 (0 O + ^kfP “ V 3 ^ 2A+e0 Q ) J ^ 
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L A (6,i) = 


L A (7,i) = 


L A (8,i) = 


1 r 9? . fl 

~ 2 h T Sin8 0 
S 3 n i 2 

- T ( ir> 

1 V . n / - V 2 

i vb I sin0 o ( st: } 

i 


n-Cosine Lead-Lag Equation 


C 1 L 
nc 

oc' 1 ' 11 + 

C 1 L (2,1) + 6 1 L (3,1) + B 1 L (4,1) 
ns nc 5 nc nc 7 ns nc 

+ 4 

l nc (5 - l > 

+ ^ 
ns 

L (6,1) + rj 
nc nc 

• . 

L (7,1) + rj L (8,1) 
nc ns nc 

+ 6 1 
nc 

L nc (9 - U 

+ B 1 

ns 

• 

L (10, i) + (J) 1 

L (11, i) + C 1 L (12, 
nc nc nc 

+ 6 1 
nc 

W 13 >« 

+ 0 

X 

L (14, i) + 0 
nc y 

L (15, i) +0 L (16,1) 
nc * y nc * 

+ R 

ys 

L „c <17 ’« 

+ R 
xs 

ft 

V* 18 ’ 17 + ir 

L „c (19>i) + 




^3 

L (21,1) = 0 
nc 


(3.13) 


L nc (1 ’« ' 

- \ + 

(“l - 

-2. 

“p) 

• 2 ft j 

sm 0 q H 


l 2 


2 £ 3 

-4 

r 


2 

e + n 

3 

- v T 

L nc (2 ’ 1) ' 

-n {2 

v Cd0 

a 

£ 4 

4 

+ i 3 

+ v T 

L K (3,i) = 

" \o 

{ \ 

- (5J 

l -2. 

, - V ! 


,-2 - 2 . 


6 P 0 O 


’SL 


f-2 - 2 n 


L (4 , i) 
nc ’ 


Vb -y- sin0Q n 2 

“< 2 T <6 k 0 + V v t <e o + W - v t <- 2 x + 5 e o» 


L (5,1) 
nc * 


- e k 0 { \ sin20 o } + c k 0 (i \ • “f } sin0 o COs6 0 
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1? 

-Vy\ 

L nc (6>1) ' ' \ V S T S “ 6 0 " 


T - ^ >e o-i S L> | 

L nc (8,i) - - 2 n-j- jj- 

1 


7r3 


tt4 


tt3 


r £ 5 “l 

L nc (9 > i) = {2 — <^n + 3 J ~ V 4- (0 n + ^n) - V 4" (-2 A+e 0 A ) } 7T 


kO 


4 w O ' Y kO' 

ft. 


'O' J fi, 


o 2 2 


-3 

l nc (10 ’ i) ’ l vE X sln0 O 2 " T 

i 

L nc (11 "l) * 4 VS T sln9 0 ^ 

•t3 2 

'■no 012 * 1 ) ' - T < TT > 

1 

L nc ( 13 ,i) = | V b 4j- sin6 0 (^) 

1 

L .c <14 ’« ‘ ' { n 4 <B p + B kO> + S n K 2 £ 1 ( £ >' 

1 

L nc (15 -« ‘ { n4 4 kO K 2 <4 )2 

1 

j3 

L nc (16>i) = { VT (0 O + \o } ’ 6 n T 2 XV} JT7 

x 

x 2 Q^. 2 

W 17 - 1 ’ - - 5 „t <s^> 

l n c< 18 -« ■ f ,T «k0 ( a^> 

L nc (l9 -« - - S nT \ «Fi> <^> 2 

t oc <2 °,i) - - .4 Ct n &> n, „ tt„,> R ( 4 > 2 


n 2 s kO “2 2 ,x Fi 


£ 2 - 2 

L „c< 21 '« ■ 5 „T h 2l 3 <V 'C 7 

1 


ft. 

x 
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n-Sine Lead-Lag Equation 


C 1 L (1 , i) + C 1 L (2,i) + B 1 L (3 , i) + 8 1 L (4 , i) 
ns ns nc ns ’ ns ns ' ' nc ns ’ ' 


+ (f) 1 L (5, i) + (J) 1 L (6, i) + r; X L (7,i) + r 1 L (8,i) 
ns ns T -nr» nc 7 ^Vi o -n c> 5 ' t-. v 7 


nc ns 


ns ns 


nc ns 


+ C L ns (9 ’ i) + ^ + 4>i L_(ll,i) + cl L_(12,i) 


nc ns 


ns ns 


ns ns 


+ L (13, i) + 0 L (14, i) +6 L (15, i) +0 L (16, i) 

ns ns y ns v ««= ~ 9 ' 


x ns 

ii 


x ns 

l 


+ K, L „s< 18 .« + T ’'ns* 19 ' 1 ) + T L n S <20 ' 1 > 


xs ns 


f ?3 I‘ ns ( 2 1 * i) = 0 


(3.14) 


L (1, i) 
ns * 


W 2 ’ 15 


-2 ~2 -2 2 -2 -2 

- + (cy - Wp) sin 0 Q + 4> k0 (oy - ajp) sin0 Q cos0 Q 

l 2 - I 3 2 I 4 

- — e + 3- n - v — B p e 0 

C d0 

n{2v _0 A_ + v l_ xo 0 + i SL } 


L „» (3 - i > 


“2 “2 “2 2 —2 -2 
- <f> k Q { - (<jy - Wp) sin 0 Q } - (cy - co F ) sin0 Q cos0 ( 


1 r . Q 2 

- I v “ -J- sm0 o n 

Ks «,±) - „ (6 0 + * k0 ) + v|- (-2 X + ie 0 > } 


W 5 >« 


L „s <6 - 1 > 


L (7, i) 
ns * 


" e k0 { “L " (W L * Sln2e 0 } + W“l “ W F ) sin0 O COS0 O 
£ 3 

- Vy A 

1 - 1? 

2 V b -y sin0Q n 

* “ 2 a 4 V 3 X0 O g SL } ft. 

x 
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p3 

L ns (8,l) - V 2 - (T 

1 

■ns' 9 ’ 13 ■ ( 2 T (8 p + B kO> - v T (6 o + \o> - V T ( - 2x+ ie o )} ST 


1 _ Q 

L ns (10 ’ i) = - 2 V b T Sin9 0 2 n $r 

i 


n 


1 


= - }vb A-sine 0 

l 3 °1 2 

L ns (129i) " - T^ 

1 

L ns (13 * 1 > * l vE f sl ” e 0 ( 57> 2 


i 3 


t2 


ft 


1 s2 


L ns (14 ’ i) = { 5 n ^ + B vn ) + 6„ V E, > ( 7 =r > 


n 3 


kO 


n 2 


ft. 

l 


_ ft] 2 

L ns (15,l) 6 n ~2~ C kO h 2 ( ft^ ) 

L ns (169i ) - { - 6 „ V T (6 0 + W + 6 n f 2 *>> ^ 


7:3 


n. 


ns 


ns 


n. 


^2 2 


n 2 

s li 

n 2 


T f 1 O J \ £ 

1J ns v * J ‘ ,/ “ u r 


£ 2 


sT> 

1 

ft. 

^kO ( 

ft. ; 

l 


f n 

kO n l 

VX Fi 


*1 \ 2 


ft. 


L ns (2 °- i)= - { „I fi 2 *2,x <*H> R( ^> 2 

L ns (21,l) = " 6 n T C k0 S 2 n 3 ^Fi* ( ft7 ) 

X 

Collective Torsion Equation 

*£ T o (1 -« + 6 m T c < 2 ' 1) + 4 T c <3 * 1) 

+ *M T c (4 - i) + 8 m I c (5 > 1) + T c <6,1) 

+ T c (7,i) + gj T c (8,1) + i* T c (9,1) 
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(3.15) 


+ e T (10, i) + -f- T (11, i) + 6 T (12, i) + 

y c K c y u 


T (1 , i) = T1 


T (2,1) = T2 
c 


T (3,1) = T3 
c 


T (4,i) = T4 


TT. 


T (5,1) = T5 


T c ( 6,i> - T6 ^ 

K 

I (7,1) ■ T7 ( jj- ) 


T c (S,i) - T8 ( ) 2 

T c (9,i) - T9 <^> 2 


T c (10,i) - - t F1 T18 


( 2 . 


T (ll,i) = n 2 a Fi ) T18 


5 

n. 


T 9 U2,i) - - I Fi T21 ( jji- ) 2 
I c (13,i) -n 2 U Fi > T21 <^-) 2 


Alternating Torsion Equation (For even N only) 
^ V'.i> + 4 T a ( 2,1) + 4 I A (3,i) 
4 T A (4>i> + 4 "a”’ 1 ' + ‘-M V 6 ’ 1 * 


T c (13,i) = 0 
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+-M T A <7 - 1) + 6-m I a < 8 - 1) + V 9 >i> - 0 


•M A 


T A (l,i) = T1 


T A (2,i) = T2 


T A (3,i) = T3 


T A (4 * i} " T4 sT 
1 


t a < 5 - 1) • 15 a- 


I A (6 - i) ■ 16 57 

1 

^1 2 

V 7 ’^ - T7 ( IT } 

X 

Q 1 2 

V 8 ’ 1 * = T8 ( sT> 

i 

^1 2 

t a (9 -« - T9 < 57 > 

i 


n-Cosine Torsion Equation 


C T nc (1 *« + C T nc (2 >« + £c T „c (3 -« + C T nc«*« 


+ <4 T „ c (5ll) + T -- (6 ' I) + T - (7 ' l)+ *-- T ~< 3 .« 


ns nc 


+ 8 nc T nc (9 ' 1)+ 4 + 4 T„(U,i)+ci T„(12,l) 


ns nc 


nc nc 


ns nc 


+ < "nc T nc a3 -« + 4 T nc <14 - i) + 4 V” 5 ' 1 * + 6 » V (16 >« 


+ 8 T (17, i) +0 T (18, i) +0 T (19, i) + R T (20, i) 
y nc x nc y nc xs nc * 


+ R T (21 , i) + R T (22, i) +R T (23, i) 
ys nc xs nc ys nc 
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+ T T nc (24 ’ U + T V* 25 ' 1 ’ + V 26 ’ 1 ’ 

5i t 2 

+ T T nc (27 ’ 1) + T V* 2 *' 11 + 5 3 T „c (2, -° 


= 0 


(3.17) 


T (1,1) = T1 - n T7 
nc 


T (2,1) = n T4 
nc 


T (3, i) = T2 - n T8 
nc 


T (4,1) = n T5 
nc ’ 


T (5,i) = T3 - n T9 
nc * 


T nc <M> - 


n T6 


y 2 -‘> - 


T4 ^ 


n. 


T (8,1) = 
nc 


T (9,1) = 
nc 


T nc a0 - 1> 


2nT7 57 

1 

T5 37 
1 

n i 

■ 2 ” 18 57 
1 


T (11,1) 
nc * 


= T6 


Q 1 

Cl. 

X 




T nc a3 -« 


T (14,1) 
nc * 


T (15,1) 
nc 


“ 2 n T9 TT 

^1 2 

- T7 ( IT > 

X 

n i 2 

- T8 < 57 > 

< 2 

- T9 ( 5J7 ) 

1 
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_ '“'I 

T nc (16,i) ={ 6 n T13 (1) - h 2 6 n T20(l) } ( ) 

T (17 , i) -{ 6 T14 (1) + h 6 T19(l) > ( ) 

nc n z n U . 

x 

T nc W.±) =U n T10(l) -h 2 « n 117(1) ) gi 

X 

^c* 19 ’ 1 * - t{ „ I11 < 1 > +S 2 5 n T16(1) 1 ^ 

^1 2 

T (20,i) = T19(l) 5 ( 7^ ) 

nc n 

^1 2 

T nc (21,l) T20(l) « n (jji) 2 

^1 

T (22, i) = T16(l) 6 ( 7=p- ) 

nc n 

T (23,i) = T17 (1) S ( ) 

nc n U. 

x 

T nc (24, 1 ) - ni « F1 ) 6 n T20(l) ( jji ) 2 


T „c (25 .» ‘ - h 2 12.x (t Fi> T19<1) «„ E < ^ > 
T nc (26,D - - h 2 n 3 «„) 120(1) S„ ( ^ ) 2 

X 

T nc< 27 ’» - n, a Fi ) 217(1) 6 n jx 

1 

T nc (28,i) - - h 2 n 2jX (*„> 116(1) « n R ^ 

T nc (29,l) - - h 2 n 3 « F1 ) 117(1) 6 n g 
n— Sine Torsion Equation 


2 

2 


2 


< |> 1 T (l,i) + (f. 1 T (2,i) + 3 1 T (3, i) + B 1 T (4,i) 
T ns ns * Y nc ns * ns ns * nc ns 


+ C 1 T (5,i) + 5 1 T (6,i) + (J) 1 T (7,i) + (J) 1 T (8,i) 
ns ns * ^nc ns * Y ns ns * T nc ns ’ 
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• . • . • . • , 

+ S 1 T (9,i) + B 1 T (10, i) + C 1 T (11, i) + T (12, i) 
ns ns nr ns 9 r»c n c 9 ' nc v 5 ' 


nc ns 


ns ns 


nc ns 


•• . M . , •• 

+ Cf) 1 T (13, i) + B 1 T (14, i) + c T (15, i) + 0 T (16, i) 
ns ns * ns ns ns ns ’ y ns ’ 


ns ns 


+ 6 T (17, i) + 0 T (18, i) + 0 T (19,1) + R T (20, i) 
x ns ’ y ns * x ns 9 ys ns * 


+ R T (21, i) + R T (22, i) + R T (23,1) 
xs ns ’ ys ns xs ns 9 ' 

•• M 

+ X W 24 ' 17 + x T „.< 25 .i) + 5, T „ = (26.i) 


R ns 


3 ns 


+ X T „s (27 ' 17 + X I n S <28 - 1) + h T „ s (29 . 1 > 


= 0 


(3.18) 


T ( 1 , i) = T1 - n T7 
ns ’ 


T (2,i) = -n T4 
ns 


T (3, i) = T2 -n T8 
ns 


T (4,1) = - n T5 
ns ’ 


T (5,1) = T3 - n T9 
ns 


T^, 1 ) - - n 16 

'.s* 7 ' 1 ) ■ T4 sT 


Si, 


T ns (M) - -2n T7 Jj- 

i 

Si 

W 9 - 17 ■ 15 sT 

1 

T n.( 10 « M ‘ - 2n 18 C 

X 

T (11, i) =16^ 
ns fi. 
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a (if 

T M < 12 .« - - 2 » T9 rr t 

^1 2 

T ns (13,i) - T7 ( js- ) 

T (14, i) = T8 ( ) 2 

ns 

^1 2 

T (15, i) = T9 ( o~ ) 
ns 9 

^1 2 

T (16, i) = {6 T14 (2) + h 6 T19(2) } ( rr~ ) 

ns n z n 

__ 2 

T (17, i) = {6 T13(2) - h„ T20(2) } ( tt~ ) 

ns n z n 

a. 

V 18 ’ 1 ) ' tS n T11(2) + S 2 6 n T16<2)1 

T (19, i) = {6 T10(2) - h 6 Tl7(2)} ^r- 

ns n z n 

^1 2 

T (20, i) = 6 T20(2) ( ^ ) 

ns n 

^1 2 

T (21, i) = 6 T19 (2) ( ) 

ns n M. 

ft 

T ns i22 ^ = 5 n T17(2) ft7 

T ns (23 > i} = 6 n Tl6(2) ^ 

- W T20<2) 6 n ( C )2 

V 25 ^ = - s 2 n 2 jX a Fi > 119(2) s n r ( ^ ) 2 

°1 2 

I ns ( 2 6,i) ■ - S 2 n 3 « F1 ) T20(2) « n C > 

^1 

V 27 - 1 * - \ W F1 ) niQ) «„ < 5 ^ > 

^1 

^(^.D - - h 2 n 2>x ai Fi ) T16C2) * n ^ 


41 



^s' 29 - 1 ' - - K 2 "3 ( V I17(2) 

The various coefficients for the torsional equations are: 

r2 


11 - - 5 I1 - 4 


2 ? k 0 2 I cos6 0 + 2 B p X I sln0 o 


+ ( 1 - 5 ) [ 


MB 3 
mR 2 


mR 


MB 2 , , . 2 a 2 a . 

~ 2 ~ J (sxn 0 Q - cos 0 Q ) 


+ V? k0 ( _ 4 + 3 2 e ^ + v \ 0 (1 + ? k0 )[ 4 ? k0 (P kO + *V 


r v- 

+ 3 - X + 3 - X 5 ] 


l 3 l 2 

+ vx a ( T + X 2 e > 


T2 “n ? k0 3 '’kO 


- + XT cos20 o + sin2 V - v T<o 


mR 


mR 


+ vB k0 (1 + ^ n )[ — (0 n + - 2 XC,.J - V X T cos0 


’kO 


’kO v 0 


kO 


a k 0 


0 


- V(1 + S 2 0 ){ - -7“ Ct- *t 2 S-x Cw, + s_xe„ + 


C d0 , £ 4 


I 3 


’kO kO 


0 1:0 


£ 3 


3 A (6 0 + ^kO ~ 2 ? kO 0 p " 2 ? ko\o ) 

+ ^ X ( X - 5 <9 0 + * k0 » } - v X A £ C t0 

T3 ■ ’"n e k0 - V T + T 5) • 6 k0 T ■ T *k 0 *i cose o ' T *i sin6 o 


+ (1 - i)(e p + <W t cos2e 0 + sta2e o 1 W B p + B k 0 > 


mR 


mR 


i 4 


£ 3 


S, 2 - 


+ v[ - (e 0 + \o - W P - ^kofW + T<- x+ 2 i 0 o + 2 i \c> - T i X1 

-L _3 

+ vB ko (1 + 4’ 1 X < 6 k0 + B p > % + W - 2 T x <6 P + V 1 

v2B k0^k0 ‘ ~ ^ X + "3~ 2 e) ■ T ^k0 * B k0 + B p 3 + ^'k0 J 

I 3 

3 X (0Q + 4 > k0 - 2 C k 0 B p - 2 ^ k 0 \o^ 
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T4 = 


T5 = 


T6 = 


T7 = 
T8 = 


2 —3 

+ X x (X - 5 (9 0 + »kO» > - (6 P + 8 k0> 

2 ?kO (1 - 5 > 1 si " 2e o + cos2e 0 1 + V f ? kO E 


mR 


mR 


-3 -2 

+ I v£ ? koT cose o +v \ T 5 

-2 -3 

+ vB k0 (1 + 4o ) C T X * + i s V (sin0 o + \o cos0 o )] 

+ -vb(x A ~ y) y- - y v ^ J ( T" + ®' e )~ § ST 

2 T" 0 kO + ^kO^ ^ *1 sin0 O + 0 kO^ 

- 2 (1-i) [ ^ sin 2 B 0 + ^ cos 2 9 0 ] - VC kQ ( . £ + l 4 i > 


mR 


mR 


+ Vg k0 (1 + ^kO^ ^ 4 (0 O + \o 5 kO 0 p ^kO ^kO^ " 4 ? k0 (0 p + 0 k(P 

3 —3 _2 

+ f (- 2 X + i (0 O + <t> k0 ))] -VX A ( f + y- i ) 


tt 2 - 


+ 2 ^ (1-i) [ sin 2 9 0 + 

0 . .. r Imb3 2 a . 

- 2 \o (1_e) [ TT cos e o + 


mR 


. sin0 

O 

vT( 

0 

1 

*P 

l 2 Xj 

I MB2 

2* 

] + 

I 4 

2 

mR 

COS 0Q 

V> ~4~ 

I MB2 

• 2 q 

] - 


mR 2 

sin 0 Q 

) 

v x 


- 2 (1-i) ( 


■ Sib 3 


mR 


yp- ) cos0 Q sin0 Q 


mR 

c 


’kO T T k0 ; 


’kO 


+ v8 k0 c 1 + 4o> 1 2 -T- X + T (e " + x 1 


0 kO 


3 -2 

+ PX A [ 2 ^-( e 0 + + k0 ) 


- U . S ) , I « B 3 + W 


mR 


mR 


2 m* eJ £ 1 


-3 “2 -2 

- «k 0 V - V *i “ se o + T *ko *1 sine 


0 


43 


T9 = 


T10(l) 


T10(2) 


Tll(l) 


Tll(2) 


T13(l) 


t n r r MB 3 2 Q J 

b (1-e) ? kQ [ 2 cos 0 Q + 

mR 

2 V 5 3 k0 (1 + ? k0* 1" 


“3 ~2 —2 

T 3 ko + F ^ko x i cos6 o + V *i sin0 o 


X MB2 . 2 Q 
0 sm 0 

nR 2 0 

1 -|vi 

k r 

l 3 

S k0 

3 

<sin6 0 + *10 

cos0 Q ) - 

|vb 

(X i 


COS0, 


b, r 


*■ (l-e) e,.„ [ cos 2 0„ + sin 2 0 o ] + \ V b 0,^(1 + ? 2 J ^ 0„sin0 


kO L p 2 
mR 


0 2 
mR 


J k0' 


’k(T 3 (T 


J 3 _ l 2 

2 ~ ? k0 2e ~2 ^k0 2 ~2 X I cos0 O 2 A e X i COS0 O 


I 2 » 


0 


+ 0 2< ^kO X I Sin0 O + 2 ^k0 ^ 1-e ^ t ^2 COS 6 0 + „2 Sin 9 0 ^ 


0 


*k0 


mR 


mR 


+ V3 k0 ^ + ^k0^ ^ 4 ^po®r\ “ o ^ s ^ n0 n 1 


4 *k(T0 


'0 


I 3 2 l 2 - 
= 2 T ? k0 + 2 ~2 C k0 x i cos6 0 

, o /•. r\ r X MB3 2 a , I MB2 , 2 n , U 4 - l 3 

+ 2 (1-e) [ - 5 - cos 0 Q + — sin 0 Q ] - VC k0 — - V X A y 


mR 

+ v8 k0 (1 + ? ko> [ T (e o + W - 2 (x - 5 V 1 


mR 


-3 


J? 2 

2 T ? k0 2 ~2 ? k0 x i cos0 O 

- 2 (1 -> t ^F cos20 o + ^F sin20 O J + ^kO I +vi Al 


mR 


mR 


- v \o (1 + c ko> 'T ( V \o> - ! T (1 - is o> 1 

_2 2 

= - 2 y ? k0 - 2 e -J" ^kO " 2 ~2 X I COS0 O " 2 * e *1 cos0 O 


+ T 2 ^ko x i sin0 o + 2 c ko (1_S) 1 “IT cos20 o + “IT sin20 o ] 

mR mR 

+ vB ko (1 + <£o> 1 X 9 o ? k0 - I B T sin0 o 1 


£ 3 2 


-3 -2 -2 

= " T ? ko " 3 ko (3 P + 3 ko ) V " T ? ko x i cos0 o - V (3 P + 3 ko )x i sin0 o 


44 


T13(2) 
T14 (1) 

T14(2) = 

T16(l) = 
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_2 

T21 =-?k 0 X ■ 1 *1 cos0 o + 1 *ko *i sin0 o • 1 v b e k 0 (1 + c k 0 } X sin0 o 


Rigid Body Equations 
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(3.19) 


where W W F1 + w f2 + W yN + + W g 
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It should be noted that equations (3.7) - (3.20) are valid for both 
articulated and hingeless rotors. 

For articulated blades, the following quantities should be set equal to 


F » L 

while for hingeless blades 


w r = 0 , co T = 0 , oj t1 - 0 


“T2 ° 

Depending on the type of rotor, articulated or hingeless, the appropriate 
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substitutions indicated above, have to be made in Eqs. (3.7) - (3.20) and the 
resulting equations represent the linearized stability equations. The equations 
for pitch, roll and the supporting structure elastic modes, presented in the 
following pages, are different for the case of articulated or the case of 
hingeless rotors. This difference is primarily due to the terms representing 
the transfer of moments at the hub due to the blades. In articulated rotors, 
without lag dampers or hinge springs, flap and lead-lag moments at the hinge are 
zero, whereas in the case of hingeless rotors, these moments are nonzero. Hence, 
two sets of equations have to be given. These equations are provided below and 
they are identified as applicable to "Hingeless Rotors" and "Articulated Rotors" 
respectively. 
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Pitch 
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Elastic Mode Equations of the Supporting structure 


Symmetric Bending in X-Y plane (Horizontal) 
2 „ „ x 3 


J >1 «Fi> r2{ »ic < - ■ " B P T - v \o T-»f ? to 0 o - b E T- si ” e o > 


+ T <B „ + 2B ko> - v T (2e o + V + v T < 2X - 'V + v T x > 


x2 


Is 

+ v T (B P - B k0> e o > 

+ £< * 2 ^- 2 ^ T - 4 «o> 

i 4 

+ *lc K 

3 _2 

+ ♦is < "V'T " t sin0 o > 

—3 —2 

+ B iV 12 B k 0 - * T (e o + V + v 4 < 2X - 5 V 5 5T 


kO 

t 3 


3 2 

+ b L < vt + \o 4 + 2 1 v E 4 aine o > 


5 

ft. 


*1 -2 c d0 4 4 

+ ^ c < A - 2 v -4 V - v^-X0„> 


’kO 


0 


ft. 

i 


ft. 


+ 4 -4 20 o B p > 35 


+ *J C < - ^ V b 


r a 


x 2 


sin 0 ^ > 


IT 

i 


+ ej c < | v b 


r i 


x 2 


ft 


1 ^2 


sin0 " > ( IT ) 


o 


*’i 4 ^1 2 

+ 0 ls< 0 p X > ( ftl > 


"i 4 o 

+ - v > 4 } 




" J 2 x 2 “I 2 

+ 0 < 3 4 - + 4 (g + B ) + 2 h,l > ( 4 ) z 

x p 2 2 p kO 2 iii 

. x 3 j 3 ft 

+ 0 < 3 v 4 + 3 v 4" >4 
x p 3 kO 3 ft^ 


54 


* E 2 1 ? 7 2 ^1 

+ e „ < 2 f= + V V (e„ + - v 4r 2 X > J- 


P 2 


3 '“0 ■ r kO' 

ft. 


a. 


+ R < -21 > ( -i ) 2 

ys v ft. 

•• 1 

^1 o 

+ t < - «1 < V 21 > < Sir > 

l 


ft 


+ h < h 9 H, U M ) 2 l > ( — ) 2 } 


’3 ' “2 "3 v ~Fi' - ~ • 'ft. 

l 


+ 2 n i,x w Fi> r3 1 +ii < - v t x > 

+ ^ < 2 T % * <W - v X (e o + V - v T ( - 2X + 2i V > 

_ m ±_ ^ r ^ 

a 4 


* 3 


x 3 


+ si < -2 V 4 - - V ~ xe - g__ > 


M 


o 6 sl n. 


+ ii < - T V b ^ sin 9 _ - B Z > 


M 


0 H p g ST 


ft. 


si. 


”i I _ 0 3 ft-i o 

+ 3 m K 2 vh T sin0 o > ( fiT } 

X 


-3 

+ r 1 < _ iL. 

3 


I 2 - > K i 2 

T e > ( 5 T ) 

1 


+ 9 y l n <- v T (e o + \o> +2 V T X > > 

^2 I 3 7 2 

+ — < Ho (&„.)(- v4-(0. + *,J + 2v^-X) 

r zri J u ku z 


°1 

ft. 


n. 


> 7 ^ > 1 . 

it. 1 

1 


+ 6 x < W P Z - r *l^Fl ) W F 1 - *>1 (£ F 2 > W F 2 - VV W EN " W 


W, 


UN 


2 •• 


m sbxy ^1 ^1 ” k sbxy ^1 “ 0 


r ^?2 

i-pi 


H 1 (x) Wg (x) dx > 


( 3 . 23 ) 


Symmetric Bendihg In X-Z plane (Vertical) 
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(3.27) 


Elastic Mode Equations of the Supporting Structure 


Symmetric Bending in X-Y plane (Horizontal) 
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(3.28) 


Symmetric Bending in X-Z plane (Vertical) 
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Torsion 
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In the relatively long set of equations presented in the preceding pages, 

Eqs . (3.1) - (3.6) represent the equilibrium position or trim equations, while 
Eqs. (3.7) - (3.30) are the linearized stability equations written in the rotor- 
plane or multiblade coordinates. The trim equations are the same for both 
articulated and hingeless rotors. .For the case of hingeless rotors, the stability 
equations are given by Eqs. (3.7) - (3.25) and for the case of articulated rotors, 
the stability equations are given by Eqs. (3.7) - (3.20) and (3.26) - (3.30). 

These equations can be used to analyze the aeroelastic and aeromechanical stability 
of a twin rotor system, with a buoyant envelope (Fig. 2), in hover. 

The stability equations can be written more compactly by using a matrix re- 
presentation 

[M] {q} + [C] {q) + [K] {q} = 0 (3.31) 

where [M] , [c ] , [k] are constant coefficient matrices in which the elements 

are dependent on the equilibrium quantities, and {q} is the generalized coordinate 

vector, which can be written as 



^5 

A , 


For a four bladed rotor, the value of n, in the blade equations, is 1. The 
corresponding generalized coordinates are 
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It can be seen from the equations associated with alternating modes, Eq. (3.8), 
(3.12) and (3.16), that these equations are decoupled from other degrees of freedom. 
Similarly, the equations corresponding to the other degrees of freedom do not depend 
on the alternating modes. Hence, Eqs. (3.8), (3.12) and (3.16) can be solved inde- 
pendently, thereby reducing the size of the matrices [M] , [C] and [K] . Based on this 
property of the alternating modes, Eq. (3.31) can be split into three groups of 
equations, namely: 
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[m 2 ] 

{q 5 > + [c 2 ] 

{q 5 > + 

[k 2 ] 

{q 5 > 

= 0 

(3.33) 

[m 3 ] 

{q 6 > + [C 3 ] 

{ V + 

[k 3 ] 


= 0 

(3.34) 

The order of 

[Mj] 

the various matrices is 
, [Cj], [Kj] 25 x 25 

given 

below 



[m 2 ] 

, lc 2 ], [k 2 ] 

3x3 





[m 3 ] 

After obtain 

, [C 3 ], [K 3 ] 
ing the equilib 

3x3 

rium state, the 

three 

groups 

of Equations (3.32) 


(3.34) can be solved separately for the stability analysis. The information about 
system stability is obtained from an eigenanalysis of Equations (3.32) - (3.34). 

The various results obtained together with the physical interpretation of these 
results are presented in the next chapter. 

3.2 Equations for Single Coupled Rotor/Body Model 

It is evident from the preceding discussion that the mathematical model of a 
multiple rotor system, coupled with a supporting structure, is algebraically compli- 
cated. To develop confidence in this model it seemed prudent to use it first for 
simulating the behavior of a single rotor system coupled with a fuselage. For a 
coupled rotor/body system, in ground resonance, including the effect of the aero- 
dynamic loads, high quality experimental data has been published by Bousman in Ref. 6. 
A comparison of the results obtained from the analytical model developed in this report, 
with experimental data [Ref. 6] is a reasonable approach for validating the equations. 
In this section, the equations of motion for the multirotor model, derived in the 
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previous sections are modified to study the aeromechanical stability of a single 
rotor/body model of a helicopter in ground resonance. The modifications intro- 
duced consist of deleting a number of degrees of freedom so as to simulate the 
particular configuration tested by Bousman [Ref. 6]. A brief description of this 
test is provided below. 

Bousman [Ref. 6] has obtained excellent experimental data for the aeromech- 
anical stability of a hingeless rotor on a special gimbaled support, simulating 
body pitch and roll degrees of freedom. The rotor consisted of three blades and 
five different configurations were tested. The different configurations re- 
present different blade parameters characterized by the nonrotating natural fre- 
quencies of the blade in flap and lag, pitch-lag coupling and flap-lag coupling. 
The rotor was designed such that most of the blade flexibility was concentrated 
at the root by building in root flexures. The rotor assembly was supported on 
gimbal which had pitch and roll degrees of freedom. In this report, the analyti- 
cal results obtained are compared with the experimental results, presented by 
Bousman, for rotor configuration 1, where the designation of this configuration 
is consistent with that in Bousman ? s paper [Ref. 6]. Configuration 1 had dif- 
ferent stiffnesses in flap and lag respectively, the corresponding non-rotating 
flap frequency was 3.13 Hz and that for lead-lag was 6.70 Hz. The airfoil 
cross-section of the blade was cambered and has a zero lift angle of attack 
equal to -1.5 degrees. A substantial part of the experimental data was 
obtained for zero pitch setting, however, due to the presence of camber the rotor 
produces a small amount of thrust at this pitch setting. The rotor blades were 
rigid outboard of the flap and lag flexures which were located at a radial station 
0.105R. There was no flap-pitch or pitch-lag couplings for this configuration. 
Furthermore, the blade was very stiff in torsion. In the case of the experiments 
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conducted for pitch angles other than zero, the experimental set up was so de- 
signed as to introduce the changes in pitch angle outboard of the flexures and 
hence there was no flap-lag structural coupling for these cases. The structural 
damping in body roll was very small in comparison with that for body pitch. The 
body pitch and roll frequencies were controlled by cantilever springs on which 
the gimbal was mounted. It was stated in Ref. 6 that the body pitch spring was 
selected to provide a dimensionless body pitch frequency of about 0.12 at a nominal 
rotor speed of 720 R.P.M. and the roll spring was selected to give a dimensionless 
roll frequency of about 0.28. (The frequencies are nondimensionalized by dividing 
by rotor speed.) As indicated in a letter by Bousman to the authors the design 
objectives for the model were dimensional frequencies of 1.44 Hz in pitch and 
3.36 Hz in roll. However the actual measured frequencies were 2 Hz in pitch and 
4 Hz in roll. From the experimental results presented in Ref. 6 it is evident 
that over a wide range of fl(200~l000 R.M.) the pitch and roll frequencies are very 
close to 2 Hz and 4 Hz respectively. Hence, for the present study, the pitch 
and roll frequencies are chosen to be 2 Hz and 4 Hz. With this combination of 
frequencies, at a rotor speed of 750 R.P.M. , the lead-lag regressing mode fre- 
quency coalesces with the body roll frequency causing an aeromechanical in- 
stability . 

The degrees of freedom required to study this aeromechanical stability prob- 
lem are: the fundamental flap and lag modes for each blade and the pitch and roll 

degrees of freedom of the body. In this class of problems, it has been established 
that the collective flap and lag modes do not couple with the body motion and 
thus, these inodes are not considered. Therefore, the number of degrees of freedom 
for the aeromechanical problem are six. These consist of: cyclic flap (8^ c > f3j g ), 

cyclic lead-lag 4^ g ), body pitch (0) and body roll (<f>) . The relevant equili- 

brium and stability equations, for this problem are given in the following sections. 
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3*2*1 Static Equilibrium Equations 
Flap : 


9 _? _2 2 p “3 t 2 

e 0 {w F + (<o L - co F ) sin 0 c + - 3 - + T* 
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When there is no structural flap-lag coupling, the terras containing sin0 
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and cos must be deleted from the above equations as well as in the 
c 

stability equations given below. 


0 


0 



9 


ZL 


where 0^ 

0 

c 



is 

is 

is 


the effective angle of attack 

the collective pitch setting of the blade 

the zero lift angle of attack 


3.2.2 Stability Equations 
n-cosine Flap 
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For a three bladed rotor, the value of n, in the above equations is 1. 

The inflow ratio, X, used in the calculation of the aerodynamic loads is 
taken from [Ref. 3] 

X = 

In the last equation 0 q is the effective angle of attack of the blade. 

As indicated in Ref. 6, a cambered airfoil was used in the model rotor 
tested, thus 
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where 0^ is the collective pitch setting of the blade and 0^ is the zero lift 

angle of attack. The static equilibrium equations, Eqs. (3.35) and (3.36), are 
used to evaluate the blade equilibrium positions and Eqs. (3 .37) -(3 .42) are solved tc 
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determine the aeromechanical stability of a single rotor helicopter, in ground 
resonance. 

The procedure followed in this analysis is slightly different from the 
procedure used to analyze the stability of a vehicle in hover, presented in 
Chapter 2. During hovering flight, the vehicle equilibrium conditions have to be 
satisfied. Whereas for ground resonance problems, the equilibrium condition of 
the vehicle do not have to be satisfied. Hence, the collective pitch angle of 
the blade, 0 , is not an independent variable to be evaluated from the 
equilibrium conditions of the vehicle, and it becomes a prescribed quantity. 

The procedure for the analysis of ground resonance problem is as follows. For 
a given value of the collective pitch setting of the rotor, under the prescribed 
conditions of operation, the equilibrium deflections of the blade have to be 
evaluated from the equilibrium equations of the blade, Eqs. (3 .35) -(3 .36) . 

Then, these quantities are substituted in the linearized stability equations, 

Eqs. (3.37) - (3.42), to analyze the stability of the vehicle. 

The results of this analysis together with the results obtained for the 
stability of the complete HHLA model are presented in the next chapter. 
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4 . RESULTS 


Based on the equations presented in the previous chapter, two types of 
problems were solved. First, the aeromechanical stability of a single rotor 
helicopter, in ground resonance, is analyzed and the analytical results are 
compared with the experimental results available in the literature [Ref. 6]. 

Next a detailed stability analysis of a multirotor vehicle, representing an 
HHLA (Fig. 2), in hover, was carried out. The results are presented in two 
separate sections. 

4.1 Results for the Ground Resonance Problem and Comparison with Experimental 

Data 

In this analysis, aimed at predicting the aeromechanical stability of 
a single rotor helicopter, the behavior of the model was studied at various 
values of the rotor speed £2. The results of this aeromechanical problem are 
presented in Figs. 4-8 together with the experimental results. Also presented 
in these figures are the results obtained by Johnson [Ref. 7]. The aerodynamic 
model used by Johnson was based on a dynamic inflow model, whereas a quasisteady 
aerodynamic model is used in this report. The data used for this analysis is 
presented in Appendix C. 

The variation modal frequencies with £2 are presented in Fig. 4, together 
with the experimental data obtained in Ref. 6. The progressing flap (3^) an d 
the progressing lead-lag (£ ) frequencies increase very rapidly with £2. The 
lead-lag regressing mode (£ ) frequency evaluated from our analytical model is 
in excellent agreement with the experimental results. The body pitch (0) and 
roll (())) frequencies have slightly higher values than the experimental results. 

The damping in pitch as a function of £2 is shown in Fig. 5. The analytical 
results are in good agreement with the experimental data. The variation of 
damping in roll as a function of £2 is shown in Fig. 6. It is evident that for 
this case the analytical results yield values which are somewhat higher than 

the experimental data. The differences observed between our analytical results 
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and the experimental points, for the frequency and damping in body modes, could 
be explained as follows. In our calculations, the numerical values used for 
the stiffness and structural damping in body pitch and roll modes are evaluated 
based on pitch frequency equal to 2 Hz and roll frequency equal to 4 Hz . 

Fig. 7 presents the variation of damping in lead-lag regressing mode with 
The results of the present analysis show slightly better agreement than the 
results obtained in Ref. 7 with inflow dynamics. It is also important to note 
that in the region, beyond 800 R.P.M., our results are in excellent agreement 
with the experimental results, while the theory with inflow dynamics predicts 
higher values. 

Changes in damping of the lead-lag regressing mode as a function of the 
collective pitch setting of the blade are presented in Fig. 8. At ft » 650 R.P.M., 
the results shown in Fig. 8a indicate that the theoretical analysis used by 
Bousman [Ref. 6] predicts a much lower value for the damping than the experimental 
results. The present analysis shows considerably better agreement. At ft = 900 
R.P.M., the experimental results indicate a lead-lag regressing mode which is 
always stable, but the theoretical results shown by Bousman [Ref. 6] imply an 

instability which becomes stronger beyond a collective pitch setting of 2 degrees. 
As evident from Fig. 7b, the results of our analysis predict the correct trend 
and the predicted damping levels are much closer to the experimental results. An 
item to be noted in these figures (8a, 8b) is that the curve representing our 
analytical results starts from an angle 0^ = - 1.5 degrees. Although Fig. 8 
contains an experimental data point corresponding to 0^ = - 3 degrees, we have 
not computed the results for this pitch setting. 


c 
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The above comparison shows good agreement between our analytical results 
and the experimental results for the aeromechanical stability of a helicopter in 
ground resonance. Therefore, it can be concluded that our analytical model for 
the dynamics of the coupled rotor/vehicle system and the method of solution for 
the stability analysis are valid. 

Finally it should be noted that more comprehensive results comparing the 
experimental data [Ref. 6] with the results from the mathematical model 
developed in this report, using a quasisteady aerodynamic model, can be found 
in Ref. 8. Additional results showing the sensitivity of the results to un- 
steady aerodynamic effects was presented in Ref. 9. 

4.2 Results for Multirotor Model of an HHLA 

Based on the equations presented in Chapter 3, two computer programs were 
developed to analyze the trim and stability of the twin rotor vehicle with a 
buoyant envelope shown in Fig. 2. The results are presented in three main sections. 
The first section gives the data and certain preliminary calculations for various 
frequencies. The second section presents the results of a parametric study in 
which certain relevant physical parameters of the system are varied so as to 
determine their effect on the stability of the vehicle. This parameter variation 
was also utilized for identifying the physical meaning of the various eigenvalues 
obtained in the analysis. The last section presents the physical interpretations 
of the results. These calculations were done on a vehicle without a sling load. 
4.2.1 Data for the Multirotor Model 

The data used* for the calculation of equilibrium or trim state and stability 
of the vehicle are given on the next two pages. 
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Blade data 


The HHLA model (Fig. 2) has identical rotors. 

Type of rotor: Articulated rotor 

Number of blades N 4 

Blade chord c = 2b 41.654 cm (1.3666 ft.) 

Hinge offset e 30.48 cm (1 ft.) 

Rotor radius R 8.6868 m (28.5 ft.) 


Blade precone 


6 . 


0 


Distance between elastic center and aerodynamic center 
Distance between elastic center and mass center X 


0 

0 


Mass/unit length of the blade m 


7.9529 kg/m (0.1661 slug/ft) 


Principal mass moment of inertia of the blade/unit length 


MB 3 


MB 2 




-2 


1.1503 x 10 kg.m (2.586 x 10 slug ft) 


-,-3 


-3 


Aerodynamic data 
Blade airfoil 
Lift curve slope 
Density of air 

Blade profile drag coefficient 
Rotor R.P.M. 

Solidity ratio 
Xock number 


'd0 


ft 

a 

Y 


6.6723 x 10 kg.m (1.5 x 10 slug ft) 

NACA 0012 
2tt 

1.2256 kg/m 3 (0.2378 x 10 -2 slug/ft 3 ) 

0.01 

217.79 R.P.M. 

0.0622 

10.9 


Nonrotating blade frequency parameters (Articulated blade) 


,Flap frequency parameter 


UU = ( T ) 

F mlT 


Lead-lag frequency parameter 


h k 

\ - < -j ^ 

mR 
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Torsional frequency parameter 
(articulated) 

•ri - < % 

mR 

0 


Torsional frequency parameter 
(hingeless) 

“T2 ' < ^3 >* 

(Assumed) 1.895 

rad/ sec 

Damping in flap 

8 sf 

0 


Damping in lead-lag 

8 sl 

0 


Damping in torsion 

8 st 

0 


Vehicle data 




Weight of fuselage F 

W F1 

3.5919 

x 10 4 N(8075 lb) 

Weight of fuselage 

W F2 

3.5919 

x 10 4 N(8075 lb) 

Weight of underslung load 

W UN 

6.6723 

x 10 4 N(1.5xl0 4 lb) 

Weight of envelope 

W EN 

8.5539 

x 10 4 N(1.923xl0 4 lb) 

Weight of supporting structure 

w s 

9.4302 

x 10 3 N(2120 lb) 

Weight of passenger compartment 

V 

6.6723 

x 10 3 N(1500 lb) 

(Treated as a lumped structural 

load 



attached at the point 0 g on the 

structure (Fig. 2) 



Buoyancy on the envelope 


1.3748 

x 10 5 N(30907 lb) 

Geometric data 




Distance between origin 0 g and 

f a 

r l rl 

-21.946m (-72 ft) 

Distance between origin 0 and 

s 

F 2 ^F2 

21.946m (72 ft) 

Distance between origin 0^ and 




underslung load (Assumed) 

h l 

-15.24m 

(-50 ft) 

Distance between centerline and 




rotor hub 

h 2 

2.591m 

(8.5 ft) 

Distance between centerline and 




center of volume of envelope 


14 . 64m 

(48.03 ft) 
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Distance between center line and 


C.G. of the envelope 8.544m (28.03 ft) 

Structural Dynamic Properties of the Supporting Structure 

The supporting structure is modelled as an elastic structure with three normal 
modes of vibration: two normal modes for bending in vertical and in horizontal 

plane and one mode for torsion. The two bending modes are symmetric modes and 
the torsion is an anti-symmetric mode. It was assumed that the envelope and 
the underslung load are attached to the supporting structure at the origin 0 g . 

The data given above shows that the vehicle is symmetric about Y-Z plane. Further- 
more due to the presence of a heavy mass attached at the center (0 g ) of the support- 
ing structure, the mode shapes in bending and torsion for each half of the model 
are assumed to be the modes of a cantilever with a tip mass. 

Modal Displacement at F-^, F 2 and 0 S 

The symmetric mode shape in bending for each half of the supporting structure 
can be written as [Ref. 11, Page 140] 

Tij (^) = 6(^) 2 -4 + (Bending in X-Y plane) (4.1a) 

and 

n 2 ( L ) = 6 ( ) 2 - 4 ( Yl ) 3 + ( L )* (Bending in X-Z plane) (4.1b) 

where X is the coordinate of any section of the supporting structure from origin 

0 and L is the length of the supporting structure, L = 21.946m (72 ft), 
s 

Bending in X-Y plane 

The modal displacement at any location on the supporting structure during the 
symmetric bending in X-Y plane can be obtained from Eq. (4.1a). The modal dis- 
placement 


at location F^ 


3.0 

at location F^ 

W 

3.0 
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0 


at origin 0 


V°s> 


The slope due to the modal displacement, at any section, can be obtained by 


differentiating Eq. (4,1) with respect to (— ) . The slopes due to the modal dis- 

L 


placement 


at location F, 


at location F, 


d(i) hi 

Ju 

1 dn i 

n i,x U F 2 ) = L d( X ) £ f2 
L 


-0, 1823 


0. 1823 


at location 0 


1 '1 

n l X ( 0 s ) = L X 

* L d( f ) 0 

ti £ 


Bending in X-Z plane 

The modal displacement and slopes due to the modal displacement at any 
location on the supporting structure, during symmetric bending in X-Z plane can 
be obtained from Eq. (4.1b). 

The modal displacement 

at location F^ : ^2 ^fP 3.0 

at location F^ : 3.0 

at location 0 : r) o (0 ) 0 

s 2 s 

The slopes due to modal displacement 


at location F, 


n2 ’ xWn> ■ i d( | > hi 

Li 


-0. 1823 


at location F r 


1 '2 

n 2 x U F 2 ) = T X i 

’ F2 - L d( | ) F 2 


0. 1823 


at location 0 


1 '2 

n 2 x ( °s ) = L T~ 

,x s L d( — ) 0 

L £ 
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Torsion 


The mode shape for torsion, for each half of the supporting structure is 
[Ref . 10 , Page 99] 

n 3 ( l ) = sin ^ ( — ) (4.2) 

The modal displacement due to torsion 


at location : ^^fP -1.0 

at location F 3 : 1.0 

at location 0 : Do(0 ) 0 

s 3 s 

Generalized mass and stiffness data 

^ 1 . 

Generalized mass ( M ) and generalized stiffness ( K ) for the i mode of 
vibration of the supporting structure is defined as 


» - /“ f: 

Jo 

rl 


F2 2 J 
m r\ dx 


and 


K - 0)7 M 

l 

.th 


where co^ is the i modal frequency 
is the i^ mode shape 

and m is the mass/unit length (for bending inodes) , or m is the mass moment 
of inertia/unit length (for torsion modes) 

Bending in X-Y plane (horizontal) 


generalized mass M 

generalized stiffness K 

Bending in X-Z plane (vertical) 

generalized mass M 

generalized stiffness K 


SBXY 


SBXY 


SBXZ 


SBXZ 


6.801 x 10^ kg (4.66 x 10^ slug) 

7.96 x 10 7 kg/sec^(5.454x!0^slug/sec^) 


6.801 x 10^ kg (4.66 x 10^ slug) 

7 .96xl0 7 kg/sec^(5 .4 54x1 0^ slug/ sec^) 
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Torsion 

generalized mass 1.936 x 10^ kg.m 2 (l .428xl0^slug ft 2 ) 

generalized stiffness K 7 .202xl0 6 kg .m 2 / sec 2 (5.312xl0^slug 

, 2 / 2 ^ 
ft /sec ) 

Rotary inertia data for the vehicle 

The rotary inertia tensor of the various masses of the system are added 
together to obtain the rotary inertia of the complete vehicle about X-Y-Z axes. 

The inertia tensor is an assumed quantity 

T = 6.44 x 10 5 kg.m 2 (4.75 x 10 5 slug. ft 2 ) 

■ 2.59 x 10^ kg.m 2 (1.91 x 10^ slug. ft 2 ) 

Z YX X XY = ° 

To facilitate distinction between data which was available and data which 
had to be assumed, the list of assumed data is provided below: 

(1) Torsional frequency of the blade 

(2) Principal moments of inertia of the blade 

(3) The mode shapes of the supporting structure and hence the generalized masses 
and stiffnesses 

(4) Inertia tensor of the vehicle 
4.2.2 Preliminary Calculations 

In the preliminary calculations, the frequencies of various modes are cal- 
culated using elementary structural dynamics. These calculated frequencies are 
useful in identifying the various eigenvalues obtained in the stability analysis. 
Supporting structure frequencies 

Bending : If the supporting structure is considered to be a free-free beam (Fig. 9a) 

with uniform properties, the first elastic mode frequency in bending is [Ref. 10, 
page 80] . 
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0 ) 


up 2 ( ^_H 


where L is the length of the beam 

El is the stiffness of the beam 
m is the mass/unit length of the beam 
In the present analysis 

m = 21.908 kg/m (0.4576 slug/ft) 

EI = 1.471 x 10 10 N.m 2 (3.56 x 10 10 lb ft 2 ) 
L = 43.891 m (144 ft) 


thus 

( 1.51TT ,2 . 1.471xl0 10 M 
W ^ 43.891 ' ' 21.908 ' 


302.7 rad/ sec 


For the case of the vehicle being considered, a heavy mass is attached at 
the center of the beam. This heavy mass is due to the envelope and under slung 
weight. Therefore, the model for the supporting structure becomes a beam with a 
heavy mass in the center (Fig. 9b). It is shown in Ref. 10 that if the ratio between 
the mass fixed at the center of the beam to the mass of the beam is greater than 
3, then frequency of the beam in the symmetric modes becomes close to the natural 
frequency of a cantilever beam with length equal to half the length of the free-free 
beam. In the present case, even with the envelope mass alone, the ratio is 


8 5539xl0 4 

— : r = 9.07. Hence the first symmetric mode for one half of the structure 

9 .4302x10“* 

can be assumed to be the fundamental mode of a cantilever. The natural frequency 


of a cantilever in fundamental mode is [Ref. 10, page 77] 


0) = ( 0 >59 7 tt 2 El h 

c v L ’ v m ' 

where L = 21.946m (72 ft) 

Therefore, the natural frequency is 0)^ = 189.27 rad/sec (4.3) 

In addition to the heavy mass attached to the center, there are also two masses, 
representing the helicopter, attached at the two ends of the beam. Thus an 
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equivalent approximate model would consist of two cantilevers with tip mass 
(Fig. 9c). The natural frequency for this approximate model is calculated using 
Rayleigh's quotient. Since the model is symmetric about the center, only one half 
should be considered for the frequency evaluation in fundamental mode. 

The fundamental mode shape for a cantilever is [Ref. 9, page 140] 


X N 2 


X n 3 


X A 


T ^(t ) “ 6 ( L ) 4 ( L ) +( L ) 

The generalized mass for the fundamental mode is 


M - f 

J o 


mdxri 


+ 

g L 


w. 


FI 


(4.4) 


= 2.311 mL + 9 — 

g 

where m is the mass/unit length of the beam 

L is the length of the cantilever 

\J . is the weight of the tip mass 
F 1 

The corresponding generalized stiffness in fundamental mode is 

2 

K = 2.311 mL 

where U) is the fundamental frequency of the cantilever without tip mass, which 
c 

in the present case 189.27 rad/sec (Eq. 4,3), Thus the fundamental frequency 
of a cantilever with tip mass is 




2,311 mL 


W. 


2.311 mL + 9 


Fl 

g J 


where 


m = 21 .908 kg/m 
L = 21 .946 m 


m = 189,27 rad/ sec 
c 


3.5919 x 10 
Fl 


N 


96 



(4.6) 


Therefore w =0.180 co =34.18 rad/ sec 
B c 

This is the bending frequency of the supporting structure in both X-Y and X-Z 
planes in fundamental symmetric mode. The mode shape is given in Eq. 4.4 for 
one half of the structure. 

Torsion : 

If the supporting structure is considered as a uniform beam (Fig. 10a) 
then the fundamental torsional frequency of the beam is given by [Ref. 11, page 


193] 


TT , GJ . 

l ( T> 


h 


where L is the length of the beam 

GJ is the torsional rigidity of the beam 
I is the moment of inertia/unit length about center of twist 
In the present case 


So w = 


L = 43.891 m (144 ft) 

GJ = 6.4054 x 10 7 N m 2 (1.55 x 10 8 lb ft 2 ) 

I = 140.972 kg m (31.706 slug ft) 

7T . 6.4054 x 10 7 .It 
43T89T C 140.972 > = 48.25 rad/ sec 


(4.7) 


Because of a large mass attached at the center of the beam (due to envelope and 
underslung load), the model in Fig. 10a can be modified as shown in Fig. 10b. 

In this model, the beam is assumed to be clamped at the center. The natural 
frequency for the fundamental mode in torsion is [Ref. 10, page 99] 

7 r , GJ 

“ ‘ 2L <T> 


In this case L = 21.946 m (length of the cantilever beam) 
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Thus 


0 ) = 


IT 


2 x 21.946 


6.4054 x 10' 
140.972 


) 2 = 48.25 rad/sec 


(4.8) 


Actually, this value is the same as that obtained in Eq. (4.7) because the tor- 
sional frequency of a beam with length £ with one end fixed and the other end 
free is the same as that of a beam with length 2£ with both ends free, this is 
due to the fact that when vibrating in its fundamental mode the center of the 
free-free beam is a nodal point. 

In the vehicle model shown in Fig. 2, there are two helicopters attached 
to the end of the supporting structure. They can be idealized to two tip masses 
having rotary inertia which are attached to the beam (Fig. 10c). Due to sym- 
metry only one half of the model has to be considered when evaluating the 
natural frequency. Assuming the mode shape to be [Ref. 10, page 99] 

n(f ) = sinf ( £ ) (4.9) 

the generalized mass is 

m = f in 2 ( f - ) dx + 1 ) 

J Q L H L 

= IL/2 + I_, 


where I Fl is the inertia of the helicopter attached at the end of the beam. The 
generalized stiffness is 


K 


IL 

2 


2 

“nt 


where is the natural frequency in fundamental mode without tip mass. In 
the present case, = 48.25 rad/sec, and thus the fundamental frequency of the 
beam in torsion including the effect of tip mass is 


W T = 


IL 2 

2 ^T 

— + I 

2 FI 
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where I = 140.972 Kg m 
I 


8.135 x 10 3 Kg m 2 (6000 slug ft 2 ) 


FI 

L = 21 .946 m 


and 


oo. 


“t = 


NT 


(1 + 


2 x 8.135 x 10 J 
140.972 x 21.946 


= .3997 Wjyj, = 19.29 rad/sec 


(4.10) 


The fundamental mode shape and the corresponding natural frequency in torsion 
are given by Eq. (4.9) and (4.10). 


Assumption regarding the torsional frequency of the blade 

The following calculation shows why the nonrotating torsional frequency 
parameter of the blade is assumed to be = 1-895 rad/ sec. 

Torsional frequency of a blade with root spring K, is 


<f> 


K . 


0 ), 


I MB3 R 


where t ^ ie 111338 moment of inertia/unit length of the blade. Assuming 

U)^ = 6fl where (3 is the angular speed in R.P.M. 


K , = 36 Q 


I MB3 R 


Using the values 


and 


X MB3 = 1 * 1503 x 10_1 K & m 
R = 8 . 6868 m 


(a)t 2 becomes 
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OJ, 


T2 


K. Is 

(-|) ! 

mR 


36ft IktuoR i 

- ( 

mR 


where m = 7.9529 Kg/m 

ft = 22.807 rad/sec 

Thus 

, 36 x 22.807 x .11503 x 8.6868 

“ T2 = ( — 3 > 

iZ 7.9529 x 8.6868 J 


= 1.895 rad/sec 

This value of oo,^ provides a torsional frequency 6ft for the blade. 

Roll Frequency of the Vehicle 

The roll frequency of the vehicle is evaluated based on the simple model 
(Fig. 11) where the force due to buoyancy is assumed to act above the C.G.. From 
Fig. 11, the equation of roll motion can be written as 


Thus 


e x + (P Z h 3 - W EN V 6 - 0 


, P Z h 3 - W E» \ ^ 
“roll ' ( ! > 


XX 


Using the values 

p| = 1.3748 x 10 5 N 

W_ M = 8.5539 x 10 4 N 

EN 

h^ = 14.640 m 

h^ = 8.544 m 

I = 6.44 x 10 5 Kg m 2 

xx 


(4.11) 


“roll ‘ ( 


1,3748 x 10 5 x 14,640 - 8.5539 x 10^ x 8,544 M 

r ' 


6.4401 x 10“ 


= 1.4108 rad/sec 


(4.12) 
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Pitch Frequency of the Vehicle 

The pitch frequency of the vehicle, is evaluated using the same assumption 
used in determining the roll frequency, thus 


p z h 3 " w EN h 4 h 

“pitch ‘ ( I > 


yy 


fi o 

where I = 2.59 x 10 Kg m 

yy 

and ai , , = 0.7036 rad/sec 

pitch 


(4.13) 


4.2.3 Summary of the Various Frequencies 

For the sake of convenience, the various frequencies, needed during the analysis 

of the vehicle, are summarized below. These frequencies are nondimensionalized 

with respect to rotor speed ft, which is equal to ft = 22.807 rad/sec. 

Rotor Blade; (In uncoupled modes) 

Rotating flap frequency for an articulated blade is 

-2,3 e 
1 + 2 R-e 

where e = 0.3048 m 
R = 8.6868 m 

TVni c 


and 


Wp = 1.0545 


<u fl = 1.027 

P 


(A. 14) 


Nondimensional rotating lead-lag frequency is 


-2 _ 3 e 
\ 2 R-e 


Thus 


= 0.0545 


go =0.233 

C 


(4.15) 
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Nondimens ional rotating torsional frequency is 


sj- 

<p 


K, 


W 8 


+ 1 


Thus 


Vehicle : 


= 37 


GO = 6,08 
<P 


Rigid body translation 


oo =0 
RX 

oo =0 
RY 


(4.16) 


Rigid body rotation 

Pitch 

00 .*, = 
pitch 

0.7036 oriOC 10 -1 

22.807 -3085 x 10 

(4.17) 


Roll 

C0 - - = 

roll 

1.4108 ,, nr , rt -l 

22.807 * 6185 X 10 

(4.18) 

Supporting structure 

flexible 

modes 



Bending in X-Y plane 

^ S BXY 

34.18 

22.807 

= 1.499 

(4.19) 

Bending in X-Z plane 

^ S BXZ 

34.18 

22.807 

= 1.499 

(4.20) 

Torsion 

W ST 

19.29 

22.807 

= 0.846 

(4.21) 


4.2,4 Equilibrium (Trim) Results without Sling Loads 

An equilibrium analysis for the vehicle in hover is performed using 
the data given in the previous sections and assuming that the magnitude of 
the underslung load is zero. 


Total weight of the vehicle W = W„_ + W_ + W„ + W c > 

EN S H S 
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= 8.5539 x 10 4 + 9.4302 x 10 3 + 2 x 3.5919 x 10 4 + 6.6723 x 10 3 
= 1.7348 x 10 5 N (39000 lb) 

Buoyancy of the envelope 

= 1.3748 x 10 5 N (30907 lb) 


Weight to be supported by the two rotors 
= 0.36 x 10 5 N (8093 lb) 

Thus each rotor has to develop a thrust = 0.18 x 10"* N (4046.5 lb) 
The various equilibrium values for these conditions are: 
Equilibrium flap angle of the blade 

lead-lag angle 
torsion angle 

Collective pitch angle 

Inflow ratio 

Thrust developed by each rotor 

Thrust coefficient 


■ 2.302 degrees 
^k0 = -3.963 degrees 
\o = “0.115 degrees 
9 q = 4.206 degrees 
X = 0.03272 


= 0.1797 x 10 N (4040 lb) 


C T = 0.00158 


Buoyancy ratio 


BR = 


1.3748 x 10~ 
1.7348 x 10 


= 0.792 


As indicated previously the equilibrium values are evaluated using an 
iterative procedure. Therefore, the difference in thrust equal to 30 N is a 
very small quantity which is assumed to represent a converged value. This 
quantity will change the equilibrium angles only in 4th or 5th decimal point. 
4.2.5 Stability Results 

Using the equilibrium values from Section 4.2.4, a stability analysis 
was performed. From the stability analysis, the eigenvalues of the linearized 
system of equations are obtained. Since the present model consists of 31 
degrees of freedom, one obtains 62 eigenvalues. Before proceeding to obtain 


103 



the stability boundaries, the eigenvalues have to be identified. To 
identify these eigenvalues, a parametric study was performed, in 
which the stiffness of the supporting structure and the rotary inertia of the 
vehicle were varied. The various other input quantities were kept fixed. Nine 
cases, listed below, were studied. It should be noted that for all these cases, 
the trim quantities are the same because the trim quantities are independent of 
the parameters modified in the parametric study. 

Case 1: Data as presented in the previous section. 

Case 2: The generalized stiffness in torsion, K , is increased from 

2 ^2 

7.202 x 10 6 Kg ^ to 1.21 x 10 7 Kg m ■ . This increases the 

sec sec 

torsional frequency of the supporting structure, from 19.29 rad/sec 

to 25 rad/sec. In nondimensional form, the increased torsional f re- 
25 

quency is 0) ST = — = 1.096. 

Case 3: Torsional frequency of the supporting structure is = 1.096. Rotary 

inertia of the vehicle is increased in pitch and roll I is increased 
5 2 6 2 

from 6.44 x 10 Kg m to 2.0 x 10 Kg m . I is increased from 

yy 

2.59 x 10^ Kg to 4.7454 x 10^ Kg m^. 

Case 4: Bending stiffness of the supporting structure is increased in both dir- 

ections. The generalized stiffnesses K S gxY an< ^ K SBXZ are increase< i tL ' om 
7 2 8 2 

7.96 x 10 Kg/sec to 1.7 x 10 Kg/sec . This increases the bending 
frequency of the supporting structure in both directions from 34.18 


rad/sec to 50 rad/sec. In nondimensional form, co 


SBXY W SBXZ 22.807 


2.192. 

Case 5: Torsional frequency of the supporting structure is increased to 40 


rad/sec. The generalized stiffness K corresponding to this frequency 

o JL 

7 2 2 

is 3.098 x 10 Kg m /sec . In nondimensional form 
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W ST 

2.192. 


40 _ 

227307 = 1.754. The bending frequencies are = a> gBXZ 


5 2 6 

Case 6: Rotary inertia in roll is increased from 6.44 x 10 Kg m to 2.0 x 10 

Kg m , aj gT - 1.754, ^ SBXY = “sBXZ = 2 ’l 92 * 

Case 7: Rotary inertia in pitch and roll are increased 

I - 2.0 x 10 6 Kg m 2 , I = 4.7454 x 10 6 Kg m 2 
xx yy 

W ST = 1 * 754 ’ “SBXY = “SBXZ = 2,192 

Case 8: A spring is introduced in the X-direction of the translational motion 

such that the nondimensional X- translational frequency is 3 =0.01. 

Kx 

Also I = 2.0 x 10 6 Kg m 2 , I = 4.7454 x 10 6 Kg m 2 , S__ = 1.754, 

xx yy st * 

"sBXY’ = ^SBXZ = 2,192 

Case 9: A spring is introduced in the Y-direction of the translational motion 

such that the nondimensional Y-translational frequency is = 0.01. 

I = 2.0 x 10^ Kg m^, I = 4.7454 x 10^ Kg m^, 3 = 1.754, 

xx yy 9 ST 

“SBXY = ^SBXZ = 2 ' 192, 

The results of the stability analysis for these nine cases are presented 

in Tables I, II and III, each column representing one case. For convenience 

a row number is also used on the left hand side of the Tables I, II and III. Thus 

tin t In 

(I,j) refers to the eigenvalue in I column and J row. 

It has been previously noted that the alternating mode of the blade is 
independent of the other degrees of freedom. Thus there are two sets of identical 
eigenvalues (presented from rows 28 - 33) one for each rotor. These are 

- 0.5200 ± i 0.5845 x 10 1 

- 0.6562 ± i 0.7265 

- 0.6522 x 10' 2 ± i 0.2346 

Since it was assumed that the torsional frequency of the rotating blade was 
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gj^ = 6,08 (Eq. 4.16), the eigenvalue - 0.5200 ± i 0.5845 x 10 should correspond 
to alternative torsion mode. The second eigenvalue - 0.6562 ± i 0.7265 
corresponds to alternating flap mode. This mode has a high damping and the 
damped flap frequency is 0.7265. The other eigenvalue -0.6522 x 10 ± i 0.2346 

corresponds to the alternating lead-lag mode. This mode has a low damping. 

These three modes are all damped modes. These alternating modes have the same 
values for all the nine cases. Since the alternating modes are stable and re- 
main unchanged for all cases, no further discussion of these modes is presented. 

It can be seen from the results in column 1, there are 5 eigenvalues 
with frequencies close to 0.7 (1, 15-1, 19) of which one eigenvalue has a 
positive real part (1, 17). The eigenvalue is 0.1024 ± i 0.7428. This eigen- 
value can correspond either to the torsional frequency of the structure or low 
frequency progressive lead-lag mode or collective flap mode. Because the tor- 
sional frequency of the structure (shown in preliminary calculations) is = 

0.846 (Eq. 4.21) and the collective lead-lag frequency is 0.233 (Eq. 4.15), 
the progressive low frequency lead-lag mode could be close to the torsional 
frequency of the supporting structure. Hence there can be coupling between 
these modes. In order to identify the various eigenvalues, the parametric study 
was performed with an aim to decouple various blade and vehicle modes. 

Consider the results for case 7 (column 7 in Table II). First all the 
eigenvalues will be identified as shown on pp . 107-108 and subsequently discussions 
of each mode are given in Section 4.2.6. Since the HHLA model (Fig. 2) consists 
of two rotor systems, the stability analysis will provide a pair of eigenvalues 
for each rotor degree of freedom. 

It is easy to identify the blade torsional, flap and lead-lag frequencies. 


From the preliminary calculations, the torsional frequency is 0 )^= 6 . 08 (Eq. 4.18), 

the flap frequency is 1.02 7 (Eq. 4.14) and the lead lag frequency is U)^=0.233 

(Eq. 4.15). 


106 


Torsional Frequency of the Blade 


Collective 

-0.5198 ± i 0.5845 x 10 1 


(7,5) 


-0.5199 ± i 0.5845 x 10 1 


(7,6) 

High Frequency or 

1 



Progressing 

-0.5207 ± i 0.6846 x 10 


(7,3) 


-0.5202 ± i 0.6845 x 10 1 


(7,4) 

Low Frequency or 

i 



Regressing 

-0.5207 ± i 0.4845 x 10 


(7,7) 


-0.5202 ± i 0.4845 x 10 1 


(7,8) 

Flap Frequency of the 

Blade 



Collective 

-0.6341 ± i 0.7361 


(7,15) 


-0.6534 ± i 0.7210 


(7,16) 

High Frequency or 

1 



Progressing 

-0.6555 ± i 0.1726 x 10 


(7,9) 


-0.6558 ± i 0.1727 x 10 1 


(7,10) 

Low Frequency or 




Progressing 

-0.6565 ± i 0.2737 


(7,20) 


-0.6562 ± i 0.2737 


(7,21) 

Lead-lag Frequency of 

the Blade 



Collective 

-0.6210 x 10 -2 ± i 0.2337 


(7,22) 


-0.6524 x 10 -2 ± i 0.2346 


(7,23) 

High Frequency or 


1 


Progressing 

-0.1676 x 10 ± i 0.1136 x 

10 

(7,13) 


-0.6923 x 10 -2 ± i 0.1252 x 

10 1 

(7,14) 

Low Frequency or 

_2 



Progressing 

-0.4893 x 10 ± i 0.7536 


(7,18) 


-0.8424 x 10 -2 ± i 0.7772 


(7,19) 
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Rigid Body Modes 


Rigid body translation in X and Y directions 

0.0 0.0 (7,1) 

0.0 0.0 ( 7 , 2 ) 

0.9833 x 10~ 5 ± i 0.2333 x 10 -4 (7,27) 

Rigid body rotation 

pitch -0.3621 X 10 -1 0.0 (7,24) 

-0.1446 x 10 -1 0.0 (7,26) 

roll -0.5947 x 10 _4 ± i 0.3510 x 10 _1 (7,25) 


Elastic Modes of the Supporting Structure 
Bending in X-Y plane (Horizontal) 

-0.1029 x 10 -2 ± i 0.2173 x 10 1 (7,11) 

Bending in X-Z plane (Vertical) 

-0.6136 x 10" 2 ± i 0.2188 x 10 1 (7,12) 

Torsion -0.6372 x 10 -2 ± i 0.1782 x 10 1 (7,17) 


4.2.6 Interpretation of the Physical Meaning of the Eigenvalues 
Blade Torsion Modes 

It is assumed that the uncoupled torsional rotating natural frequency of the 

blade is co. = 6,08 (Eq. 4,16) . Thus the eigenvalues corresponding to this fre- 

<P 

quency must represent the collective torsional mode frequency. Frequencies 
corresponding to co^ ± 1 represent the cyclic mode frequencies. The cyclic 
modes are <p^ c and cj)^. The progressing mode has a higher frequency and the 
regressing mode is the lower frequency. All these modes have negative real 
part indicating a positive damping and hence these modes are stable. 
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I 


j Blade Flap Modes 

The uncoupled rotating natural frequency in flap is shown to be = 1.027 

Eq. (4.14). In the presence of aerodynamics, the flap mode is heavily damped. 

Thus the damped natural frequency in flap should be less than 1 .027. In the 

present case, the damped flap frequencies are 0.7361 and 0.7210. These fre- 

P 

quencies correspond to the collective flap modes. Frequencies corresponding 

to 0) D ± 1 are the cyclic flap modes. These modes are also heavily damped. In 
P 

this case, both cyclic modes are progressing modes, one with higher frequency 
and the other with a lower frequency. (When the collective mode frequency is 
less than 1, then both cyclic mode frequencies are progressing modes. Ref. 3.) 

Blade Lead-Lag Modes 

The uncoupled rotating natural frequency in lag is = 0.233, Eq. (4.15). 

This particular frequency will appear in the eigenvalues as a collective 
lead-lag frequency. Another typical property of lead-lag mode is that these 
modes are very lightly damped. Since the collective lead-lag frequency is less 
than 1, both the cyclic mode frequencies are progressing. All these three 
I modes are lightly damped. 

Inspection of the eigenvalues reveals that all the blade modes are 
associated with two sets of eigenvalues. This is caused by the presence of 
the two rotors each with its own set of blade modes. 

Rigid Body Translation Modes 

There are four eigenvalues corresponding to the rigid body translation 
in X and Y directions. Two of them having zero real and zero imaginary parts. The 
other eigenvalue set has a very small positive real part and a very small imaginary 


I 
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-5 -4 

part. (0.9833 x 10 ± i 0.2333 x 10 ). The reason for choosing this eigen- 

value as one corresponding to the rigid body translation is given below. 

Comparing the results given in columns 7 and 8, it can be seen that all 
the eigenvalues except a few remain the same. The results of column 8 are 
obtained by introducing a translational spring in the X-direction. The spring 
constant is prescribed to yield a natural frequency of oscillation m - 0.01. 
From the results of column 8, it can be seen that the second eigenvalue corres- 
ponds to this frequency having a value of 0.5621 x 10 ^ ± i .1172 x 10 

* 

The eigenvalues corresponding to motion must be 0.0 + i 0.0, (8,1) , and 
0.9762 x 10 -5 + i 0.0, which is assumed to be equivalent to zero. Then 
motion results in a pure translatory motion. A similar observation can be made 
when a translational spring is introduced in the Y-direction (results of 
column 9) leaving the translational motion in X-direction free. It is seen 
from the results that there is an eigenvalue corresponding to a frequency 


,-5 


.-1 


= 0.01 which is 0.3349 x 10 ± i 0.1171 x 10 (9,2). The eigenvalues 

corresponding to translational motion in the X— direction becomes 0.0 + 
i 0.0 (9,1) and 0.9784 x 10~ 5 + i 0.0 (9,27). Then R^ motion becomes a pure 


translational motion. 

The previous statements also imply that when the R^ motion is oscil- 
latory, Ry motion becomes pure translational motion and vice versa. How- 
ever, when both R and R are free (results of column 7) the combined R , R 
x y x y 

motions have eigenvalues which are complex conjugates. This oscillatory mode 

is a divergent mode, but the frequency and damping are very small. This 

indicates that R and R motions cannot be separated, 
x y 

Recall, as indicated on P. 105 (l,J) stands for the I t ^ 1 column and J t ^ 1 row. 
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Rigid Body Rotation Modes 

There are four eigenvalues corresponding to pitch and roll motions. These 

are 

,-4 . . „ .... . „-l 


-0.5947 x 10 ± i 0.3510 x 10‘ 

,-l 


-0.3621 x 10 A ± i 0.0 
,-l 


(7,25) 

0, 2b) 

-0.1446 x lO'- 1 ± i 0.0 (7,26) 

The oscillatory mode corresponds to the roll mode and the other two pure 
damped modes correspond to pitch mode. These statements are further clarified 
by discussion presented below. Using Eq. (4.11), the roll frequency of the vehicle 
is 

( P Z h 3 ~ W EN h 4 ,5s _1 
roll '‘I ' 

XX 

Substituting the various quantities 


03 


roll 


a , 1.3748 x 10 5 x 14.64 - 8.5539 x IQ 4 x 8.544 1 


2.0 x 10 


22.807 


Thus 


03 


0.8006 n „ in in -l 
= 0.3510 x 10 


roll 22.807 

This calculation shows that the eigenvalue -0.5945 x 10 ^ ± i 0.3510 x 10 ^ 

(7,25) corresponds to the roll mode. 

Using the same elementary Eq .( 4 . 1 1) y the pitch frequency is 

( 1.3748 x 1 0 5 x 14.64 - 8.5539 x IQ 4 x 8.544 M 1 
Pibch ‘ 4<y454 x 10 6 22.807 


0.5197 

22.807 


= 0.2279 x 10 


-1 


But this frequency is not evident in the eigenvalues. Note that for a tandem 
rotor system, the pitch mode is a heavily damped mode. When the damping 
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is in excess of the critical damping, the pitch motion becomes a pure damped 

motion. However when the inertia in pitch is increased, this mode will also 

become a oscillatory mode [Ref. 3]. The reason for the presence of a relatively 

high damping in pitch can be explained using Fig . 2. For positive pitching motion of 

the vehicle, rotor system R^ moves up and rotor moves down. If 0^ is the pitch rate, 

then the rotor R.. has an upward velocity of H 0 , is experienced by rotor R. . This in- 
1 r i y i 

creases the net inflow velocity sensed by rotor R^ . If the net inflow is increased, the 
effective angle of attack experienced by a typical blade sect ion decreases . This in effect 
decreases the thrust developed by the rotor R^ . Similarly for rotor system R^, 
the net inflow velocity decreases which in effect increases the angle of attack 
and hence the thrust. The combined effect of the increase in thrust for rotor 
system and decrease in thrust for due to a positive pitch rate 0 , tends 
to restore the vehicle to its equilibrium position. This restoring force depends 
on 0^ and produces damping in pitch. When this damp ing is high, the pitch motion be- 
comes a pure damped motion. In the present case, the damping in pitch is suf- 
ficiently high so that the eigenvalues have only negative real part. 

It is well known that for second order system with damping above the critical 
damping, an increase in inertia will bring the two eigenvalues closer provided that 
this increase in inertia is such that even with the increased inertia the system 
is still overdamped. A further increase inertia will make the eigenvalues to 
become complex conjugates. This effect is evident from the results by comparing 
the columns 6 and 7 . 

From column 6, the eigenvalues corresponding to pitch are 

-0.8143 x 10 -1 ± i 0.0 (6,24) 

-0.1183 x 10" 1 ± i 0.0 (6,26) 
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The rotary inertia in pitch for this case (case 6) is I = 2.59 x 10 6 Kg m 2 . 

6 2 

When the rotary inertia is increased to 4.7454 x 10 Kg m , keeping the other 
parameters the same, the eigenvalues corresponding to pitch motion become 
(column 7) 

-0.3621 x 10 _1 + i 0.0 (7,24) 

-0.1446 x 10 -1 + i 0.0 (7,26) 

This shows that the eigenvalues have approached each other. This validates the 
statement that the pitch mode, in this case, is a overdamped mode. 

Elastic Modes of the Supporting Structure 

In the present analysis, the supporting structure is modelled by three 
normal modes: two for bending and one for torsion. The two bending modes corres- 

pond to bending in X-Z plane (Vertical) and bending in X-Y plane (Horizontal) . 

The bending mode in X-Z plane has higher damping than that corresponding to the hori- 
zontal bending mode. The explanation is the same as that given for pitch motion, 
in previous section. 

The eigenvalues for 

bending in X-Z direction is -0.6136 x 10~^ ± i 0.2188 x 10^ (7,12) 

bending in X-Y direction is -0.1029 x 10 ^ ± i 0.2175 x 10 1 (7,11) 

torsion is -0.6372 x 10 -2 ± i 0.1782 x 10 1 (7,17) 

The identification of these modes is based on the frequencies assumed in obtaining 

the results presented in column 7 

W SBXY = W SBXZ = 2,192 

and 

“ST ' 1 - 754 
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4.2.7 Coupling of Various Modes 


The coupling between various blade modes and body modes are shown in 
Figs. 12-16. Since the HHLA model vehicle (Fig. 2) consists of two rotor 
systems coupled by the supporting structure, it was shown in Section 4.2.5 
that the stability analysis provides a pair of eigenvalues for each rotor degree 
of freedom. Hence for the purpose of identification, in the presentation of the 
results shown in Figs. 12-16, the rotor modes will be referred to as mode 1 
and mode 2, such as collective flap mode 1, collective flap mode 2 and high 
frequency flap mode 1 and high frequency flap mode 2, etc. 

Figure 12 illustrates the variation of the eigenvalues of blade lead-lag 
modes and the supporting structure bending modes as a result of an increase in 
the bending stiffness (K 013W ) of the supporting structure in X-Y (horizontal) 

ODAI 

plane. The bending stiffness K gBXY was increased in increments from 5.09 x 10 

O 

N/m to 1.74 x 10 N/m, such that the corresponding uncoupled nondimensional 

bending frequency in X-Y plane assume ^ t ^ ie va l ues ^ggxY = ^ ^ ’ 1*^99, 

1.754, 2.192, where the frequencies are nondimens ionalized with respect to the 

rotor speed of rotation £2, where fl = 217.79 R.P.M. The arrows in the figure 

indicate the direction along which the eigenvalues of the modes change due to an 

increase in K . The eigenvalues of the other modes, which are not shown in the 
SBXY 

figure, remain unaffected by the variation in It can be seen from Fig. 12 

that the bending mode, in X-Y plane, of the supporting structure is strongly couplei 
with the high frequency lag mode 2. The high frequency lag mode 2 which was 
initially unstable becomes more stable as Kgg X y is increased. The damping in the 
bending mode in X-Y plane decreases asympototically with an increase in frequency 
and this mode is always stable. The low frequency lead-lag mode 2 shows a slight 


114 


decrease in damping as i s increased. The eigenvalues corresponding to the 

bending mode in X-Z plane and the high frequency lag mode 1 are not affected by 
the changes in However, since these two modes have nearly equal fre- 

quencies it can be seen that the high frequency lag mode 1 is unstable* 

Figure 13 presents the variation of eigenvalues of the blade lead-lag modes 
and the supporting structure bending modes as a result of an increase in the 


bending stiffness (K^^) of the supporting structure in X-Z (vertical) plane. 

6 

The bending stiffness was increased in increments from 7.96 x 10 N/m 

8 

to 1.74 x 10 N/m and the corresponding nondimensional uncoupled bending frequency 


in X-Z plane (^ SBXZ ) assumed the values ^ SBXZ = 1.499, 1.754, 2.192. It can be 
seen from Fig. 13 that the bending mode in X-Z plane is strongly coupled with 
high frequency lag mode 1 . The high frequency lag mode 1 which was initially 
unstable becomes a stable mode as K PTJV „ Is increased from 7.96 x 10 7 N/m 

odaZ, 

^SBXZ = to 1.09 x 10 8 N/m = 1.754). But a further increase in 

8 

^SBXZ t0 x 10 N/m does not affect the eigenvalue corresponding to the 


high frequency lag mode 1, indicating that these two modes are decoupled. Damping 


in the bending mode in X-Z plane decreases drastically at tne beginning and once 
the bending mode and the high frequency lag mode 1 are decoupled, the decrease 


in damping of the bending mode in X-Z plane is very small. Damping in the torsion 
mode of the supporting structure and low frequency lag mode 1 are slightly affected 


aS *SsBXZ i ncrease d* Since the torsion mode and the low frequency lag mode 1 
have frequencies which are close to each other, the figure clearly indicates that 
the lag mode 1 is unstable. The eigenvalues corresponding to the rest of the 
modes are unaffected. 


115 



Figure 14 shows the eigenvalue variation in the rotor lead-lag modes 
and the torsion mode of the supporting structure as a result of an increase 
in the torsional stiffness (K gT ) of the supporting structure. The torsional 
stiffness, K gx , was increased in increments from K gT = 1.59 x 10 6 N.m to 
3.99 x 10 7 N.m and the corresponding uncoupled nondimensional torsional fre- 
quency (w gT ) of the supporting structure are w gx = 0.4, 0.55, 0.846, 1.096, 1.2, 
1.3, 1.4, 1.5, 1.754, 2.0. It is evident from the figure that the low frequency 
lag mode 2 and high frequency lag mode 2 remain unaffected during the variations 
in K gT and these modes qre stable. In Fig. 14, the different curves are divided 
into three segments represented by points A, B, C and D. The curves between 
points A to B refer to the range of K gx = 3.01 x 10^ N.m to 7.20 x 10 N.m 

(d> = 0.55 to 0.846); the curves between points B to C refer to the range 

ST 

K_,_ = 7.20 x 10 6 N.m to 1.685 x 10 7 N.m (u> CT = 0.846 to 1.3); and th rves 
ST 

between points C to D refer to the range K gx = 1.685 x 10 N.m to 3.1 x 10 N.m 
(w gT = 1.3 to 1.754) . 

It is evident from Fig. 14 that in the ranee A to B, as the torsion^ 
stiffness K CT is increased, the torsion mode of the supporting structure becomes 
increasingly stable and its frequency is increasing; the low f reque^ zy lag mode 1 
becomes increasingly unstable with its frequency slightly increased* This clearl 
indicates that the torsion mode is strongly coupled with the low frequency lag 
mode 1* The high frequency lag mode 1 experiences a slight increase in frequency 
but its damping remains almost the same. In this range, A to B, the eigenvalues 
of these three modes have been distinctly identified based on their uncoupled 
nondimensional frequencies. In the range B to C, as the torsional stiffness 
is increased, the damping in the low frequency lag mode 1 decreases and its fre- 
quency tends to increase towards 1.0. At the same time, the damping in torsiona. 
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mode of the supporting structure decreases drastically and a slight change in 
the frequency is observed (i.e. the frequency initially increases and then 
decreases) . The high frequency lag mode 1 shows an increase in frequency with 
no appreciable change in damping. In this range B to C, the eigenvalues of 
these three modes do not exhibit a direct one to one correspondence to the 
uncoupled nondimensional frequencies, implying that all these modes are coupled. 
Hence in this range, B to C, the reference to the various modes, as torsion mode, 
low frequency lag mode 1 and high frequency lag mode 1, is only for the convenience 
of explaining the variation of the eigenvalues. When the torsional stiffness 
K was increased still further, i.e. the range C to D, the eigenvalues start 
exhibiting a correspondence to the nondimensional uncoupled frequencies indicating 
that these three modes are slowly getting decoupled. In this range, C to D, the 
torsional mode of the supporting structure has low damping and it tends to decrease 
asymptotically while the frequency increases from 1.5 to 1.75. The high frequency 
lag mode 1 shows an increase in the frequency and the mode becomes stable at 
the point D. The damping in the low frequency lag mode 1 decreases while the 
frequency undergoes a slight reduction. Beyond the point D i.e. for 
K ST ~ 3.1 X ^ t * ie ei g enva l ues °f l° w frequency lag mode 1 and high frequency 

lag mode 1 show negligible change and the damping in torsion mode remains the same 
but its frequency increases. Beyond point D all the three modes are stable. 

Another interesting observation which can be made from Fig. 14 is associated 
with the effect due to the increase in torsional stiffness K gT . When K gT is in- 
creased from 1.685 x 10 7 N-m to 3.99 x 10 7 N.m (curve in the range C to D and 
beyond) , the eigenvalues corresponding to the high frequency lag mode 1 tend to 
approach the eigenvalue corresponding to the high frequency lag mode 2 (which 
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remains unaffected during the variation in K^) and similarly the low fre- 
quency lag mode 1 approach to the low frequency lag mode 2. This behavior 
seems to indicate that as the torsional stiffness of the supporting structure 
is increased the coupling between the rwo rotors due to the torsional deformation 
of the supporting structure is eliminated. As a result of this lack of coupling, 
the eigenvalues corresponding to the high frequency lag modes 1 and 2 and low 
frequency lag modes 1 and 2 approach each other. It should be noted that elimin- 
ation of the coupling of the rwo rotors, due to the torsional deformation of the 
supporting structure, does not imply that the two rotors are totally decoupled. 

The rotors are still coupled through the bending deformation of the supporting 
structure and rigid body pitch motion of the vehicle. The presence of this coupling 
causes the eigenvalues of the low frequency and high frequency lag modes to approach 
each other rather than coalescing. 

The last observation which can be made using Fig. 14 is that the high frequency 
lag mode 1, low frequency lag mode 1 and torsion mode of the supporting structure 
undergo a reversal in their characteristics as K is increased from 1.59 x 10^ 

O i- 

N.m. Thus, the mode which was initially a distinct torsion mode becomes a low 

frequency lag mode 1; the low frequency lag mode 1 becomes a high frequency lag 

mode 1 and the high frequency lag mode 1 becomes a torsion mode. For low and 

6 — 

high values of the torsional stiffness (i.e., K ^ < 1.59 x 10 N.m (Wg T - 0.4) 
and £ 3.10 x 10^ N.m (u)^ > 1.754)) the torsional mode of the supporting 
structure, the low frequency by mode 1 and high frequency by mode 1 are all 
stable. For intermediate values of the torsional stiffness of the supporting 
one of the lag modes is unstable. 
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The variation of the eigenvalues of the collective flap modes and body 
pitch mode due to increase in body inertia in pitch is presented in Fig. 15. It 
is evident from the figure that the pitch mode is a pure damped mode. An 
increase in pitch inertia causes the eigenvalues, corresponding to the pitch 
mode, to approach each other. The eigenvalues of the collective flap mode 2 
tend to approach the eigenvalue of the collective flap mode 1. The pure damped 
nature of the pitch mode is associated with the presence of two rotors. During 
pitch motion the net inflow in the two rotor system changes. If in one rotor 
system the net inflow increases, then in the other one the inflow decreases and 
vice versa. These changes in inflow results in changes in the thrust in the two 
rotor systems. The rotor system which moves up, during pitch motion, experiences 
a reduction in thrust due to the increased inflow and the rotor system which 
moves down produces more thrust due to the decreased inflow. These changes in 
the thrust tend to restore the vehicle to its equilibrium position. Since this 
restoring force is proportional to the pitch rate, this mechanism produces a 
damping in pitch. In the present case, the pitch motion is overdamped. Hence 
an increase in inertia causes the eigenvalues, corresponding to the pitch mode 
to approach each other, as shown in Fig. 15. 

Figure 16 illustrates the variation of eigenvalues corresponding to the low 
frequency lag mode 2 and body roll mode as a result of an increase in inertia 
in roll. An increase in roll inertia tends to decrease the damping in roll, 
furthermore its frequency is also reduced. The low frequency lag mode 2 tends 
to become more stable. The roll mode, for the model vehicle, is a damped 
oscillatory mode. This is different from the pure damped mode normally observed 
in a conventional tandem rotor helicopter. The reason for this oscillatory nature 
of the roll mode is due to the presence of the buoyancy of the envelope. 
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For all the cases analyzed, it was found that the flap and torsional modes 
of the rotor are always stable. The eigenvalues corresponding to the cyclic 
flap modes and all the torsion modes are not affected by the variation in the 
quantities used in this parametric study. The alternating modes of the rotor 
were also found to be stable. The degree of coupling, as well as the relative 
strength of the coupling, between the various blade modes and the body modes 
is presented in a qualitative manner in Table IV. It is evident from this 
table that the supporting structure elastic modes are strongly coupled with 
the low frequency and high frequency lead-lag modes. 

4.2.8 Effects of Buoyancy on the Stability of the Vehicle 

The effects of varying the buoyancy ratio on the stability of the vehicle 
were also studied, by performing the stability analysis at different buoyancy 
ratios. During this analysis, only the buoyancy ratio was varied while the 
rest of the blade and vehicle parameters were kept fixed. The vehicle para- 
meters are the same as those used in Case 7, presented in Section 4.2.5. 

The results of these analyses for different buoyancy ratios are presented in 
Tables V and VI. 

Table V presents the results of the equilibrium (trim) analysis for various 
buoyancy ratios. It can be seen that as the buoyancy ratio is decreased, the 
thrust coefficient of the rotors (C T ) increases. The equilibrium angles of the 
blade in flap, lead-lag and torsion, the inflow ratio and the collective pitch 
angle, also increase with decrease in buoyancy ratios. Table V also presents 
the nondimen sional roll frequency of the vehicle at different buoyancy 

ratios. These roll frequencies are calculated using Eq. (4.11). These fre- 
quencies will be helpful in identifying the roll mode in the stability analysis. 
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Table VI presents the results of the stability analysis at different buoyancy 
ratios. The results of the stability analysis, presented in Table VI, are also 
shown in a graphical manner in Figs. 17 and 18. Figure 17 depicts the variation 
of eigenvalues of the supporting structure elastic modes with decrease in 
buoyancy ratio. The direction of arrows in the figure indicates the variation 
of the eigenvalues as a result of the decrease in buoyancy ratio. The frequencies 
corresponding to the supporting structure elastic modes are not affected by the 
variation in buoyancy ratio. However, the damping in bending in X-Y plane 
increases, the damping in X-Z plane decreases, while the damping in torsion 
mode increases. 

Figure 18 presents the variation of eigenvalues of pitch and roll modes with 
buoyancy ratio. As the buoyancy ratio is decreased, one of the eigenvalues 
corresponding to the pitch mode decreases while the other eigenvalues increases. 

The pitch mode remains a pure damped mode. The roll mode which was initially 
a stable mode becomes unstable for buoyancy ratios BR £ 0.6. The results shown 
in Table VI also indicate that when the buoyancy ratio is decreased, the damping 
in the lead-lag modes of the rotors increases while the damping in flap and torsion 
modes of the rotoes decreases. However changes in the buoyancy ratio have only 
a minor effect on the frequencies of the rotor modes. From the results shown in 
Table VI, it can be seen that for a 40% reduction in buoyancy ratio, the damping 
in torsion modes decreases by 12%; the damping in flap modes decreases by 12% 
and the damping in lag modes increases by 200%. 

The rigid body translation mode is stabilized as the buoyancy ratio is 
decreased. 

The most important observation from these results is that for buoyancy ratio 
BR = 0.7, all the eigenvalues have negative real parts indicating that the 
vehicle is stable at this buoyancy ratio. 
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5 . CONCLUDING REMARKS 


This report presents the equilibrium (trim) equations and linearized stability 
equations for the dynamics of the coupled rotor/vehicle system in hovering flight. 
The stability equations are written in multiblade (or rotor plane) coordinate 
system. Two types of problems are solved. First, the aeromechanical stability 
of a helicopter in ground resonance is analyzed, and the analytical results are 
compared with the experimental results available in literature. It was found 
that the results of the present analysis compare very well with the experimental 
results. This indicates that the theoretical model for the coupled rotor/body 
dynamics appears to be accurate. 

Next, the aeromechanical stability of an HHLA type vehicle in hover was 
analyzed. The vehicle consisted of two rotors, a buoyant envelope and an under- 
slung load attached to a flexible supporting structure. For this vehicle, the total 
number of degrees of freedom is 31 and there are 31 coupled equations representing 
the dynamics of the system. Two computer programs were developed to analyze 
the trim and stability of the vehicle. The restuls of a sample problem are also 
presented in this report. 

Before describing the conclusions obtained from the stability analysis of 
the HHLA type vehicle conducted in this study, it is important to emphasize 
that the vehicle model used in this study has only two rotors and not four rotors, 
which are present in the HHLA type vehicle under construction. Furthermore, no 
lead lag dampers were included in the treatment of the blade lead-lag dynamics 
Incorporation of such dampers would have probably stabilized any instability 
observed in the lead-lag degrees of freedom of the vehicle. 
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The stability analysis yields 62 eigenvalues, corresponding to the 31 degrees 
of freedom. The primary aim was to identify the 62 eigenvalues and relate them 
to the various modes of the rotor/vehicle assembly. This identification was 
accomplished by performing a parametric study in which the primary parameters 
allowed to vary were the bending and torsional stiffness of the supporting 
structure combined with the rotary inertia of the vehicle in pitch and roll. 

This parametric variation was done in order to decouple the blade modes from the 
vehicle and the supporting structure modes. In total, nine cases were analyzed. 

In these cases, the underslung load was not included. Based on the results 
obtained for these cases, the various eigenvalues and the coupling among different 
modes were identified and physical insight on the dynamics of the vehicle was 
developed. The most important results of this study are summarized below. 

Cyclic lead-lag modes of the rotors couple strongly with the pitch, roll 
and bending in two orthogonal plane and torsion of the flexible supporting 
structure. This shows that the frequencies of vibrations of the supporting 
structure must be separated from the frequencies of the rotor lead-lag modes. 

This also implies the importance of modeling the supporting structure with an 
adequate number of elastic modes. 

The stability analysis of the coupled rotor/vehicle dynamics illustrates 
the aeroelastic stability of the rotor, coupled rotor/vehicle aeromechanical 
stability such as air resonance and the vehicle stability in the longitudinal 
and lateral planes. Complete information of these ingredients are all captured 
by the analytical model representing the coupled rotor/vehicle dynamics. 

In the discussion of the results it has been noted that the pitch mode 
of the vehicle is a pure damped mode while the roll mode is a stable oscillatory 
mode. The oscillatory nature of the roll mode can be attributed to the presence 
of the buoyant envelope. 
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The present analysis yields a divergent mode which corresponds to a pure 
translational mode. It is also found that when the vehicle is free in 
longitudinal and lateral translational motions, the results indicate a os- 
cillatory mode for the rigid body translation. However, when the translational 
motion in one direction is restrained, the translational motion in the other 
direction becomes a pure divergent motion. This indicates that the longi- 
tudinal and lateral dynamics cannot be separated in the analysis of coupled 
rotors/vehicle dynamics for vehicles of the type considered in this study. 

The stability of the vehicle was also studied at various buoyancy ratios 
and it was found that at a particular buoyancy ratio, the eigenvalues corresponding 
to all the modes have negative real part indicating that the vehicle is stable 
at this buoyancy ratio. 

Based on the numerical studies conducted in this report it appears that the 
consistent analytical model, for the dynamics of coupled rotor/vehicle system, 
developed in this study is a valid mathematical model. The stability analysis 
of coupled rotor/vehicle dynamics yields useful information on both the aero- 
elastic stability, aeromechanical stability and also the vehicle dynamic stability. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


Case 1 

cjj *=0) =1.499.(0 =0.846 

SBXY SBXZ * ST * 

I - 6.44 x 10 5 Kg m 2 

I ■= 2.59 x 10 6 Kg m 2 

yy 

Case 2 

w sbxy =w sbxz =1 ' 499,w st =1 * 096 

I = 6.44 x 10 3 Kg m 2 

XX ° 

I = 2.59 x 10 6 Kg m 2 
yy 6 

Case 3 

a) SBXY =a) SBXZ 3l ‘ 499,a) ST =1 * 096 
I = 2.0 x 10^ Kg m 2 

I >= 4.7454 X 10 6 Kg m 2 

yy h 

Case 4 

0) SBXY =a) SBXZ =2 * 1 92 ’ W ST = 1 * 096 
I xx = 6.44 x 10 5 Kg m 2 
I yy = 2.59 x 10 6 Kg nt 2 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

-.5206 ± i .6846 x XO 1 

-.5206 ± i .6846 x 10 1 

-.5206 ± i .6846 x 10 1 

-.5206 ± 1 .6846 x 10 1 

-.5203 ± i .6845 x 10 1 

-.5203 ± i .6845 x 10 l 

-.5202 ± i .6845 x 10 1 

-.5203 ± i .6845 x 10 1 

-.5198 ± i .5845 x 10 1 

-.5198 ± 1 .5845 x 10 1 

-.5198 ± 1 .5845 x 10 1 


-.5198 ± i .5844 x 10 1 

-.5198 ± i .5844 x 10 1 

-.5199 ± i .5845 x 10 1 


-.5207 ± i .4845 x 10 1 

-.5207 ± 1 .4845 x 10 1 

-.5207 ± i .4845 x 10 1 

-.5207 ± i .4845 x 10 1 

-.5203 ± i .4845 x 10 1 

-.5203 ± i .4845 x 10 1 

-.5202 ± i .4845 x 10 1 

-.5203 ± i .4845 x 10 1 

-.6556 ± i .1727 x 10 1 

-.6556 ± 1 .1727 x 10 1 

-.6556 ± i .1727 x 10 1 

-.6555 ± i .1727 x 10 1 

-.6556 ± i .1727 x 10 1 

-.6556 ± i .1727 x 10 1 

-.6556 ± i .1727 x 10 1 

-.6558 + i .1727 x 10 1 

-.4722xlO -2 ±i. 1555x10* 

-.4722xlO~ 2 ±i. 1555x10* 

~.4693xl0~ 2 ±i. 1555x10* 

1030xl0~ 2 + i.21 75x10* 

-. 89 15xlO -2 ±i. 1474x10 1 

,2525xl0'*±i. 1467x10* 

.2525xl0"*±i. 1467x10* 

-. 6269x1 0~ 2 +i. 2189x10* 

-.1395xl0 -1 ±i. 1408x10* 

-.4803xl0"*±i. 1462x10* 

-.4803xl0~*±i. 1462x10* 

-.1143xl0 _1 +l. 1424x10* 

-. 302 lxlO -2 ±i. 11 95x10 1 

-.3021xl0~ 2 ±i. 1195x10* 

-.2948xl0 -2 ±l, 1 193x10 1 

-.7002xl0 _2 ±l. 1254x10* 

-.6157 ± 1 .7452 

-.6157 ± i .7452 

-.6341 ± i .7361 

-.6157 ± 1 .7452 

-.6483 ± i .7167 

-.6483 ± 1 .7167 

-.6483 ± 1 .7167 

-.6534 ± 1 .7210 

.1024 ± i .7428 

.7266 x 10"* ± i .8399 

.7266 x 10" 1 ± 1 .8399 

.7352 x 10“* + i .8396 

-.5095 x 10 _2 i i .7555 

-.5095 x 10 _2 ± 1 .7555 

-.5375 x 10" 2 ± 1 .7566 

-.4631 x 10 _2 ± 1 .7525 

-.1057 ± i .7306 

-.7611 x 10~* ± i .8151 

-.7611 x 10"* ± i .8151 

-.7701 x 10”* ± i .8150 

-.6565 ± i .2739 

-.6565 ± i .2739 

-.6565 ± 1 .2737 

-.6565 ± 1 .2739 

-.6562 ± i .2739 

-.6562 ± i .2738 

-.6562 ± 1 .2738 

-.6562 ± 1 .2738 

-.5697 x 10~ 2 ± i .2331 

-.5697 x 10~ 2 ± 1 .2331 

-.6210 x 10" 2 ± 1 .2337 

-.5697 x lO” 2 ± i .2331 

-.6527 x 10” 2 ± i .2346 

-.6527 x 10” 2 ± i .2346 

-.6527 x 10" 2 ± i .2346 

-.6524 x 10” 2 + i .2346 

-.8143 x 10" 1 0.0 

-.8143 .x. lO* 1 0.0 

-.3621 x 10”* 0.0 

-.8143 x 10”* 0.0 

-.2106xl0~ 3 ±i. 6174x10"* 

-.2106xl0~ 3 ±i. 6147x10"* 

-. 5949x1 0~ 4 ±i. 35 10x10”* 

-. 2 105xl0” 3 ±i. 6174x10”* 

-.1183 x 10" 1 0.0 

-.1183 x 10"* 0.0 

-.1446 x 10 * 0,0 

-.1183 x 10~* 0.0 

.9833xl0" 5 ±i.2333xl0" 4 

.9833xlO' 5 ±i.2333xlO“ 4 

. 9833xlO” 5 ±i . 2333xlO” 4 

.9833xl0” 5 ±i.2333xl0” 4 

-.5200 ± 1 .5845 x 10 1 

-.5200 ± 1 .5845 x 10 1 

-.5200 ± 1 .5845 x 10 1 

-.5200 ± i .5845 x 10 1 

-.6562 ± i .7265 

-.6562 ± i .7265 

-.6562 ± 1 .7265 

-.6562 ± i .7265 

-.6522 x 10~ 2 ± i .2346 

-.6522 x 10~ 2 ± i .2346 

-.6522 x 10” 2 ± i .2346 i 

-.6522 x 10" 2 ± i .2346 

-.5200 ± 1 . 5845 x 10 1 

-.5200 ± i .5845 x 10 1 

-.5200 ± i .5845 x 10 1 

-.5200 ± 1 .5845 x 10 1 

-.6562 ± i .7265 

-.6562 ± 1 .7265 

-.6562 ± i .7265 

-.6562 ± 1 .7265 

-.6522 x 10“ 2 ± i .2346 

-.6522 x 10 -2 ± 1 .2346 

-.6522 x 10~ 2 ± 1 .2346 

-.6522 x 10 -2 i 1 .2346 


Table I Results of Stability Analysis for Various Configuration 
Parameters 

BR = 0,792, C T = 0.00158 
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Case 5 

^SBXY^SBXZ* 2 • 1 9 2 » w sx * 1 • 7 54 
I - 6.44 x 10 5 Kg m 2 

XX 

I - 2.59 x 10 6 Kg m 2 

yy 

Case 6 

“sBXY'^SBXZ" 2 • 192 ,l °ST = 1-754 
I - 2.0 x 10 6 Kg m 2 

XX ° 

I * 2.59 x 10^ Kg m 2 

yy 

Case 7 

“sbxy““sbxz- 2-192 -“st =1-754 

I = 2.0 x 10 6 Kg m 2 

XX ° 

I ■ 4.7454 x 10^ Kg 

yy 

Case 8 

“sbxy-“sbxz- 2- * 92 '“st“* -754 - 

“rx‘ 0-01 

I - 2.0 x 10 6 Kg m 2 

XX ° 

I » 4.7454 x 10^ Kg m 2 

yy 

1 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

2 

0.0 0.0 

0.0 0.0 

0.0 0.0 

.5621x10 -3 ±i. 1172xlO _1 

3 

-.5207 ± i .6846 x 10 1 

-.5207 ± i .6846 x 10 1 

-.5207 ± i .6846 x 10 1 

-.5207 ± i .6846 x 10 1 

4 

-.5203 ± i .6845 x 10 1 

-.5202 ± i .6845 x 10 1 

-.5202 ± i .6845 x 10 1 

-.5202 ± i .6845 x 10 l 

5 

-.5198 t i .5845 x 10 1 

-.5198 ± i .5845 x 10 1 

-.5198 ± i .5845 x 10 1 

-.5198 ± i .5845 x 10 1 

6 

-.5198 ± i .5844 x 10 1 

-.5198 ± i .5844 x 10 1 

-.5199 ± i .5845 x 10 1 

-.5199 ± i .5845 x 10 1 

7 

-.5207 ± i .4845 x 10 1 

-.5207 ± i .4845 x 10 1 

-.5207 ± i .4845 x 10* 

-.5207 ± i .4845 x 10 1 

8 

-.5203 ± i .4845 x 10 1 

-.5202 ± i .4845 x 10 1 

-.5202 ± i .4845 x 10* 

-.5202 ± i .4845 x 10 1 

9 

-.6555 t i .1726 x 10 1 

-.6555 ± i .1726 x 10 1 

-.6555 ± i .1726 x 10 1 

-.6555 ± i .1726 x 10 l 

10 

-.6558 t i .1727 x 10 1 

-.6558 ± i .1727 x 10 1 

-.6558 ± i .1727 x 10 l 

-.6558 ± i .1727 x 10 1 

11 

-. 1030xl0*" 2 ± i. 2175x10 1 

-.1032xl0 _2 ±i. 2175x10* 

-.1029xl0” 2 ±i. 2175x10* 

-.1029xl0" 2 ±i.2175xl0 1 

12 

-.6136xl0“ 2 ±i. 2188x10* 

-.6136xlO _2 ii. 2188x10* 

-.6136xl0 _2 ii. 2188x10* 

-.6136xl0” 2 ±i.2188xl0 1 

13 

-.1676xl0~ 3 ±i. 1136x10* 

-.1676xl0' 3 ±i. 1136x10* 

-.1676xl0 _3 ±i. 1136x10* 

-.1676xl0’ 3 ±i.ll36xl0 1 

14 

-.7002xl0~ 2 ±i. 1254x10* 

-6951xlO -2 ± i . 12 52x10 1 

-.6923xlO _2 ±i. 1252x10* 

-.6923xl0* 2 ±i. 1252x10* 

15 

-.6157 t i .7452 

-.6157 ± i .7452 

-.6341 ± i .7361 

-.6341 ± i .7361 

16 

-.6534 t i .7210 

-.6534 ± i .7210 

-.6534 ± i .7210 

-.6534 ± i .7210 

17 

-.6372xlO _2 ii. 1782x10* 

-.6372xlO~ 2 ±i. 1782x10* 

-.6372xl0 _2 ±i. 1782x10* 

-.6372xl0 _2 ±i. 1 782x1 O 1 

18 

-.4631 x 10 -2 ± i . 7525 

-.4830 x 10 -2 ± i . 7534 

-.4893 x 10” 2 ± i . 7536 

-.4893 x 10' 2 t i. 7536 

19 

-.8424 x 10" 2 ± i. 7772 

-.8424 x 10' 2 ± i. 7772 

-.8424 x 10" 2 ± i . 7772 

-.8424 x 10~ 2 i i.7772 

20 

-.6565 ± i .2739 

-.6565 ± i .2737 

-.6565 ± i .2737 

-.6565 ± i .2737 

21 

-.6562 ± i .2737 

-.6562 ± i .2737 

-.6562 ± i .2737 

-.6562 i i .2737 

22 

-.5697 x 10“ 2 ± i.2331 

-.5697 x 10’ 2 ± i.2331 

-.6210 x 10” 2 ± i.2337 

-.6210 x 10 -2 ± i.2337 

23 

-.6524 x 10” 2 ± i.2346 

-.6524 x 10~ 2 i i.2346 

-.6524 x 10“ 2 ± i.2346 

-.6524 x 10“ 2 ± i.2346 

24 

-.8143 x 10~* 0.0 

-.8143 x 10"* 0.0 

-.3621 x 10" 1 0.0 

-.3621 x 10' 1 0.0 

25 

-.2l05xl0 -3 ±i.6174xl0 - * 

- . 594 7xlO~ A ±i. 3510x10”* 

-,5947xl0“ 4 ±i.3510xl0“ 1 

- . 5947xlO _A ±i.3510xlO _l 

26 

-.1183 x 10~* 0.0 

-.1183 x 10~* 0.0 

-.1446 x 10 _1 0.0 

-.1446 x lO" 1 0.0 

27 

.9833xlO" 5 ±i.2333xlO" A 

.9833xl0” 5 ±i.2333xl0 _A 

. 9833x1 0 -3 ±i. 2333x10”^ 

.9762 x 10‘ 5 0,0 

28 

-.5200 ± i .5845 x 10 1 

-5200 ± i .5845 x 10* 

-.5200 ± i .5845 x 10* 

-.5200 t i .5845 x 10* 

29 

-.6562 t 1 .7265 

-.6562 ± i .7265 

-.6562 t i .7265 

-.6562 ± i .7265 

30 

-.6522 x 10" 2 ± i.2346 

-.6522 x 10' 2 ± i.2346 

-.6522 x 10" 2 ± i.2346 

-.6522 x 10" 2 ± i.2346 

31 

-.5200 ± i .5845 x 10 1 

-.5200 ± i .5845 x 10 1 

-.5200 ± 1 .5845 x 10 1 

-.5200 ± i .5845 x 10* 

32 

-.6562 ± i .7265 

-.6562 ± i .7265 

-.6562 ± i .7265 

-.6562 ± i .7265 

33 

-.6522 x 10~ 2 ± i .2346 

-.6522 x 10“ 2 ± i .2346 

-.6522 x 10~ 2 i i .2346 

-.6522 x 10“ 2 i i .2346 


Table II Results of Stability Analysis for Various Configurations 


Parameters 

BR = 0.792, C T = 0.00158 
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ORIGINAL PAGE IS 

OF POOR QUALITY 












































































Case 9 

U SBXy"“sBXz” 2 • 1 92 ’ “st” 1 • 7 54 • 

S^-o.oi 

I xx “ 2-0 x 10 K 8 m 2 
I » 4.7454 x 10 6 Kg m 2 

1 

0.0 0.0 

2 

.3349xlO _5 ±i.ll71xlO _1 

3 

-.5207 ± i .6846 x 10 1 

4 

-.5202 ± 1 .6845 x 10 1 

5 

-.5198 ± i .5845 x 10 1 

6 

-.5198 i i .5845 x 10 1 

7 

-.5207 ± i .4845 x 10 1 

8 

-.5202 + i .4845 x 10 1 

9 

-.6555 t i .1726 x 10 1 

10 

-.6558 ± i .1727 x 10 1 

11 

-.1029xl0 _2 ±i. 2175x10* 

12 

-.6i36xio -2 ±i. 2188 x 10 * 

13 

-.1676xl0 _3 ±i. 1136x10* 

14 

-.6923xl0 -2 ±i. 1252x10* 

15 

-.6341 ± i .7361 

16 

-.6534 i i .7210 

17 

-.6372xl0 _2 ±i. 1782x10* 

18 

-.4893xl0 _2 ±i.7536 

19 

-.8424 x 10 -2 t i.7772 

20 

-.6565 * i .2737 

21 

-.6562 ± i .2737 j 

22 

-.6210 x 10“ 2 ± i .2337 

23 

-.6524 x 10“ 2 ± i. 2346 

24 

-.3621 x 10“ L 0.0 

25 

-.5947xl0 _4 ±i.3510xl0 _1 

26 

-.1446 x 10' 1 0.0 

27 

.9784 x 10 -5 0.0 

28 

-.5200 ± i .5845 x 10 1 

29 

-.6562 ± i .7265 

30 

-.6522 x 10 -2 ± i .2346 

31 

-.5200 ± i 5845 x 10 _1 

32 

-.6522 i i .7265 

33 

-.6522 x 10 -2 t i .2346 


Table III Results of Stability Analysis for Various 
Configuration Parameters 
BR = 0.792, C T = 0.00158 
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MODES 

Lead-lag Modes 

Flap Modes 

Hi 

fi 

■eq. 

collec- 

tive 

freq. 

Low 

freq. 

High 

freq. 

collec- 

tive 

freq. 

Low 
freq . 

1 

2 

i 

2 

1 

2 

1 

2 

i 

2 

1 

2 

Supporting structure 
symmetric bending in 
x-y (horizontal) plane 


XXX 




XX 







Supporting structure 
symmetric bending in 
x-z (vertical) plane 

XXX 


X 


XX 




X 




Supporting structure 
torsion (antisymmetric) 

XXX 




XXX 








Body pitch 


X 


X 


X 




X 



Body roll 


X 




XX 








Legend: XXX = Strongly coupled, XX = Moderately coupled, X = Weakly coupled 


Table IV Coupling Between Various Body Modes and Blade Modes 
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Buoyancy 


Ratio 

BR 

o 

CD 

^k0 

4 k0 

^k0 

X 

C T 

0) i 

roll 

0.792 

4.206° 

2.302° 

-3.963° 

-0.115° 

0.03272 

.00158 

.3510x10' 

0.7 

5.243° 

3.209° 

-5.074° 

-0.161° 

0.03820 

.00228 

.3173x10' 

0.6 

6.259° 

4.179° 

-6.453° 

-0.236° 

0.04313 

.00304 

.2761x10' 

0.5 

7.207° 

5.142° 

-7.994° 

-0.352° 

0.04743 

.00380 

.2276x10 


Table V Equilibrium Values at Different Buoyancy Ratios 

^ST =1 ' 754, ^SBXY = “SBXZ =2 * 192 ,I yy = ‘ 4, 7454x10 kg ‘ m 

I =2.0xl0^kg.m^ 
xx e 
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CRIGiNAL RAGS tS 

OF POOR QUALITY 




W ST * 1 * 754, ^SBXY * ^SBXZ = 2 ‘ 192 
I yy = 4 • 74 54 x 1q6 Kg m 2 , I xx - 2.0 x 10 6 Kg m 2 


Blade 

Torsion 

BR * 0.792 

0.7 

0.6 

0.5 


-.4954 ± i .6861 x 10 1 






-.4572 ± i .6940 x 10 1 




-.4560 ± i .5939 x 10 1 

-.5199 ± i .5845 x 10 1 



-.4559 i i .5939 x 10 1 

-.5202 t i .4845 x 10 l 


H5BIPPP3SBRBI 


-.5207 t i .4845 x 10 1 

-.4956 ± i .4858 x 10 1 



Blade 

Flap 


-6411 * i .1721 x 10 1 




-.6402 ± i .1721 x 10 1 

-.6143 ± i .1725 x 10 1 


-.6341 ± i .7361 

-.6192 ± i .7313 

-.5925 ± i .7362 

-.5514 ± i .7547 

-.6534 i i .7210 

-.6393 ± i .7158 

-.6136 ± i .7196 

-.5736 t i .7364 

-.6565 ± i .2737 

-.6426 t i .2788 

-.6162 ± i .2750 

-.5773 ± i .2577 


-.6420 ± i .2789 

-.6172 t i .2748 


Blade 

Lead-Lag 


-.1006xl0“ 1 ±i.l252xl0 1 

-.1428xl0‘ 1 ±i.l252xl0 1 

-.1957xlO~ 1 ±i. 1253x10* 

-.1676xl0‘ 3 ±i.ll36xl0 1 


-.3490xl0” 2 ±i.ll36xl0 1 

-.6190xl0~ 2 ±i. 1 136X10 1 

-.6210 x 10~ 2 ± i.2337 

IffffffffTHWlfffTI 

-.1363 x 10” 1 ± i.2335 

-.1905 x 10* 1 ± i . 2328 



-.1389 x 10 _1 ± i . 2359 

-.1920 x 10' 1 t i. 2363 

-.8424 x 10~ 2 ± i. 7772 

-.1211 x 10' 1 ± i. 7769 


-.2313 x 10” 1 + i . 7772 

-.4893 x 10 -2 ± i.7536 

-.7608 x 10~ 2 ± i . 7529 

-.1136 x 10 _1 ± i. 7521 

-.1617 x 10" 1 ± i. 7514 

Rigid Body 
Translat ion 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

.9833xlO" 5 ±i .2334xl0~ 4 

-.5716xl0~ 5 ±i.5015xl0" 4 

-. 6929xl0 -4 + i. 1 058x1 0~ 3 

-3 -3 

-.2814x10 +1 237?vlO 

Pitch 

-.3621 x 10 _1 0.0 

-.4181 x 10” 1 0.0 

HSIi|SHg9 

-.5012 x 10 -1 0.0 

- . 1446 x 10 1 0.0 

-.1025 x 10' 1 0.0 


-.4048 x io -2 n n 

Roll 

-.5947xl0~ 4 ±i.3510xl0 -1 

-.3602xl0~ 4 ±i.3174xl0~ 1 

.1155x10 3 + i . 2765x1(1 1 

. 806Qv 1 0 -4 + i 

Supporting 
S true ture 

Bending in X-Y 
Bending in X-Z 

Torsion 

-. 1 02 9x1 0~ 2 ±i. 2175x10* 

wbhibib^hi 



- . 6136xl0 -2 ± i . 2 188x1 0 1 





HB9IRISBIRRIIRRI 



Alternating 

Modes 

Torsion 
Flap 
Lead -Lag 
Torsion 
Flap 
Lead-Lag 

-.5200 ± i .5845 x 10 1 

-.4941 ± i .5859 x 10 1 

-.4706 ± i .5889 x 10 l 

-.4557 ± i .5940 x 10* 

-.6562 i i .7265 

-.6421 ± i .7213 

-.6162 ± i .7254 

-.5759 ± i .7426 

-.6522 x 10" 2 t i . 2346 

-.9658 x 10" 2 t i. 2352 

-.1389 x 10 ^ i i . 23 58 

-.1921 x 10 * ± i.2363 

-.5200 ± i .5845 x 10 1 

-.4941 ± i .5859 x 10 1 


-.4557 ± i . 5940 x 10 1 

-.6562 ± i .7265 

-.6421 ± i .7213 


-.5759 t i .7426 

-.6522 x 10“ 2 ± i. 2346 

-.9658 x 10" 2 ± i. 2352 

-.1389 x 10 1 ± i .2358 

-.1921 x 10" 1 ± i .2363 


Table VI Results of Stability Analysis at Different Buoyancy Ratios 
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Figure 1 Hybrid Heavy Lift Airship - Approximate Configuration 



Figu 
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Figure 3 Equivalent Spring Restrained Blade Model 
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OUR ANALYTICAL RESULTS 

THEORY WITH INFLOW DYNAMICS (Ref. 7) 

□ EXPERIMENT (Ref. 6) 



n, R.P.M. 


Figure 5 Body Pitch Mode Damping as a Function of 
0 c = 0 (Configuration 1) 
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C'UR ANALYTICAL RESULTS 

THEORY WITH INFLOW DYNAMICS (Ref. 7 ) 

<C> EXPERIMENT (Ref. 6) 



n, r.p.m. 


Figure 6 Body Roll Mode Damping as a Function of ft, 0 
(Configuration 1) 
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Figure 7 Regressing Lag Mode Damping as a Function of L , 
(Configuration 1) 


sec cr , sec' 


OUR ANALYTICAL RESULTS 

THEORY (Ref. 6 ) 

O EXPERIMENT (Ref. 6) 



(a) n = 650 R.P.M. 



(b) fl = 900 R.P.M. 


Figure 8 Lead-Lag Regressing Mode Damping as a 
Function of 0 C at (a) 650 R.P.M. and 
(b) 900 R.P.M. (Configuration 1) 


21.946 m 
(c) 

Figure 9 Idealization of Supporting Structure for Bending 
Type Deformations 

(a) Free-Free Beam, (b) Free-Free Beam with 
Heavy Mass at the Center and (c) Free-Free Beam 
with Masses at the Center and at the Tips 










(c) EACH HALF 


Figure 10 Idealization of Supporting Structure for Torsion 
Type Deformations 

(a) Free-Free Beam, (b) Cantilevered Beam and 
(c) Cantilevered Beam with Inertia at the Tip 
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in Supporting Structure Bending Stiffness in X-Y Plane 
(Horizontal) 



k sbxz = 7 96 x 10? N/m ~ 1 - 74 x 108 N/m 
c3 S bxz = 1499 ~ 2192 



□ SUPPORTING STRUCTURE TORSION 

O SUPPORTING STRUCTURE BENDING 
IN X-Z PLANE (VERTICAL) 

A SUPPORTING STRUCTURE BENDING 
IN X-Y PLANE (HORIZONTAL) 

■ HIGH FREQUENCY LEAD-LAG 1 

• HIGH FREQUENCY LEAD-LAG 2 
A LOW FREQUENCY LEAD-LAG 1 

• LOW FREQUENCY LEAD-LAG 2 

Figure 13 Variation of Nondimens ional Eigenvalues of Blade 
Lead-Lag Modes and Supporting Structure Bending 
Modes with Increase in Supporting Structure Bending 
Stiffness in X-Z Plane (Vertical) 


144 



■ COLLECTIVE LEAD-LAG 1, 2 

• SUPPORTING STRUCTURE TORSION 
O LOW FREQUENCY LEAD-LAG 1 

• LOW FREQUENCY LEAD-LAG 2 
A HIGH FREQUENCY LEAD-LAG 1 

• HIGH FREQUENCY LEAD-LAG 2 


Kgj = 1.59 x 10® N.m ~ 3.99 x 10 7 N.m 
c^st * 0.4 ~ 2.0 



(T 


A - B K SJ = 3.01 x 10 6 - 7.20 x 10 6 N.m 
B - C K st = 7.20 x 10 6 - 1.685 x 10 7 N.m 
C - D K st = 1.685 x 10 7 ~ 3.10 x 10 7 N.m 


Figure 14 Variation of Nondimensional Eigenvalues of Blade 
Lead-Lag Modes and Supporting Structure Torsion 
Mode with Increase in Supporting Structure Stiff- 
ness in Torsion 
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A COLLECTIVE FLAP 2 
ACOLLECTIVE FLAP 1 
•OB0DY PITCH 


l yy = 2.59 x 10 6 kg.m 2 ~ 4.75 x 10 6 kg.m 2 



Figure 15 Variation of Nondimensional Eigenvalues of 

Collective Flap Modes and Body Pitch Mode with 
Increase in Body Inertia in Pitch 
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• BODY ROLL 

o LOW FREQUENCY LEAD-LAG 2 


l xx = 6.44 x 10 5 kg.m 2 ~ 2.0 x 10 6 kg.m 2 



Figure 16 Variation of Nondimensional Eigenvalues of Low 
Frequency Lead-Lag Mode and Body Roll Mode with 
Increase in Body Inertia in Roll 
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• BENDING IN X-Y PLANE (HORIZONTAL) 
o BENDING IN X-Z PLANE (VERTICAL) 

O TORSION 


Figure 17 Variation of Nondimensional Eigenvalues of the 

Supporting Structure Elastic Modes with Decrease 
in Buoyancy Ratio 


(a) PITCH 


J(0 



gure 18 Variation of Nondimensional Eigenvalues in (a) Pitch 
and (b) Roll Modes with Decrease in Buoyancy Ratio 
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APPENDIX A 


Transformation to Multiblade Coordinates 

For an N-bladed rotor with blades evenly spaced around, the azimuth angle 
for the blade, at any instant, can be written as 

ik = i|> + 2 tt § K = 1, N (A.l) 

r k T N 

where = fit, the nondimen sional time variable. 

Let a, be a generalized coordinate associated with any degree of freedom of the k blade, 

K. 

flap or lead-lag or torsion. Since this 0^ is associated with the blade which is 
rotating, it is called a rotating coordinate. If there are N blades, the behavior 
of all the blade in that particular degree of freedom can be represented by N 
rotating coordinates . . . 0^. By suitably choosing a transformation, these N 
rotating coordinates can be transformed to another set of N coordinates, each of 
which is associated with a specific variation of all the c^'s (rotating coordinates) 
when combined , as viewed from a nonrotating frame. This type of transformation is 
called the multiblade coordinate transformation. Basically, this transformation 
transforms the rotating coordinates into a nonrotating frame. Usually, the physi- 
cal explanation about this transformation is given only with reference to flap 
motion of the blade [Ref. 3]. 

The transformation from the rotating to the nonrotating coordinate is ob- 
tained from the following operations 

1 N 

“m “ N k=l “k 

1 N k 

Om = — kJj (-1) K (for even N only) 
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N 


“nc ‘ I kll cos(n, V °k 


N 


“ns ‘ I kh sln(n V \ 


(A. 2) 


where n = 1 , . . . .L and L = 


N-l 


for odd N 


T N " 2 C 

L = — for even N 


The inverse transformation is 

L 


°k = “m + n=l (a nc COS n \ + a ns Sin + (_1) a -M (A ' 3) 


Last term will appear only even N. The proof for this transformation can be found 
in Ref . 3 . 

This transformation, given in Eq. (A.l), looks like a truncated Fourier 

series, except for the last term. The major difference between this transformation 

and the usual Fourier transformation is that here the coefficients ol,, a , a , 

M nc ns 

a are all functions of time, whereas in the Fourier series the coefficients are 
— M 

constants. That is why sometimes these multiblade coordinates are also referred to 


as Fourier coordinates. 


DitfcrGntiatiug Eq . (A>2) wxth respec 

1 N 

“m = N k=l \ 


f « f Tttm ill 

lv l. xiuv vjy 


^ 0 t ^ o r> run o f t n f ^ 

U \lHl _L- »J U 


N 


a -M N k=l ^ 


N 


N 


(A, 4) 


a + n a ■ — , E, cos(n \p. ) a. 

nc ns N k=l r k' k 


a - n a = — . E. sin(nik ) a, 

ns nc N k=l r k' k 


Differentiating again with respect to ip 
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- ,1, sin niii, (...) n-sine operator (A. 6) 

N k=l k 

These four operators are applied to each equation representing the blade degree 
of freedom and the blade degrees of freedom are replaced by the multiblade 
coordinates using Equations (A. 2), (A. A) and (A. 5). The resulting equations 
will have the multiblade coordinates as the generalized coordinates. These 
equations represent the dynamics of the rotor as a whole as viewed from a non- 
rotating frame. 
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APPENDIX B 


Application of Multiblade Coordinate Transformation 
To Multi-Rotor System 


In Ref . 1, the blade loads are derived, in a general form for a typical 
rotor with a moving hub and the rotor is assumed to operate at a specified 
constant value of ft. When using these expressions for the calculation of blade 
loads, for two different rotor systems operating at different values of ft, a 
number of special provisions described below have to be introduced. 

First note that in the general expressions for the blade loads, various time 
derivative terms are nondimen sionalized with respect to a nondimen sional time 
ijj = ftt, where ft is the angular speed of the rotor. In a multirotor system, if 
the rotors operate at different values of ft, the nondimen sional time tJj is dif- 
ferent for different rotors, which leads to inconsistency. This problem can be 
resolved using the angular speed of one rotor, say ft^ of the rotor , as the 
reference ft. Then all the time derivative terms, that appear in the blade load 


expressions for different rotors, can be suitably modified such that the non- 
dimensional time is the same for the complete set of equations. The nondimens ional 
time will be ^ = ft^t. 

The second important item is encountered while applying the multiblade 
coordinate transformation to a multirotor system. In the n=cosine and n-sine 
transformation, the following operators are applied to the blade equations: 


N 


2^ _ cos nik ( . . .) 

N k=l 


n-cosine 


2 N 

- k Z x sin n* k (...) 


n-sxne 


(B.l) 


where = iJj + 2 — 
and ^ = ftt 

In a multirotor system, these operators are different for the blade equations in 
different rotors. In general, for the i rotor, these can be written as 
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cos n4^ (. . .) 


n-cosine 


2 

N 


N 

iii 


2 

N 



sin n4^ 


i 

k 


( 


• • • 


) 


n-sine 


(B. 2) 


where 4^5 = 4^. + 2tt — ■ 
k l N 

4> = fi.t 

i l 

Note that 4^ contains the angular speed 0, . of the i rotor. These operators 
have to be applied to the blade equations in the i th rotor system to transform 
the blade equations to the nonrotating frame. When applying Eqs. (B.2) to the time 
derivative terms for the blade degrees of freedom, the transformation given in 
Eq. (A. 4) and (A. 5) should not be used directly because the time derivative is 
taken with respect to a reference nondimen sional time, say 4^ • 

The derivation given below shows the modifications which have to be incor- 
porated, in the blade equations, both for nondimens ionalizat ion and for the 
application multiblade coordinate transformation when dealing with a multirotor 
system in which each rotor is operating at a different value of Q. 

Let be the reference angular speed which represent rotor system and 

1 th 

let (3^ b e the flap degree of freedom of the k blade in the rotor system . 

Furthermore let be the angular speed of the i*"* 1 rotor system and denote by 

the flap degree of freedom of the blade in rotor system R^. The time 

derivatives of these degrees of freedom can be written in the nondimensional 

form as. 



= n 



1 d(Si t) 





i 

n 



(B.3) 
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(B. 4 ) 


K K 

dt “i d(n,t) 


aef n. d ef 

= - n ± C fr > — I 

d^ W i dip 1 


If the general load expressions, derived for typical rotor blade, are used for 
all the rotor systems in a vehicle, then the first time derivative terms in the 
load expressions for the rotor will be with respect to and that for the 
rotor R ± will be with respect to iJj 1 . If the angular speeds in the rotor systems 
are different, then the nondimen sional time will be different for the rotors* 

For consistency, the underlined expression should be used for all the first derivative 
terms appearing in the i^ rotor load expressions. Hence, the general blade load 


expressions can be used for different rotors after multiplying the first time 

^1 2 

derivative terms by ( 7— ) and the second time derivative terms by ( tt- ) . The 

fl i L i 

nondimen sional time, in this case, is ^ (Note: All the fi’s are constants). 

The following derivation shows how this nondimensionalization affects the multi- 
blade coordinate transformation. Let a 1 and a 1 be the transformed n-cosine and 

nc ns 

n-sine degrees of freedom in the nonrotating frame and let the corresponding rotat- 
ing blade degree of freedom be aj^, for rotor system Thus 


a 


nc 


a 


ns 


9 N . . 

— - cos nUL 1 Ot . 1 

N k/=l ‘ It K 


2 y • I i i 

N k=l Sin n *k “k 


2lTk 


where ty, 1 = 4 11 + „ 

k N 


and tjj 1 = J2^t , is the rotor angular speed 

Differentiating Eq. (B. 5 ) with respect to t 

*i 2^ i i i*i 

a = sin nij\ . nSi. a. + cos nik a. 

nc N k=l r k i k r k k 


(B. 5 ) 


(B. 6) 


(B. 7 ) 
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where (*) = -7 ( ) 

Equations (B.7) can be written as 


• i. i 2 i * 3 

a 1 4- nfl. a 1 = — . E cos nil ;. 1 a. 3 

nc x ns N k=l r k 1 


(B. 8 ) 


Nondimensionalizing the time derivative with respect to ip = f2jt, where £2^ 
is the reference angular speed, Eq* (B. 8 ) becomes 


* i 2 v i *i 

52- a + n£2. a = — E, cos ntk £2 a 

1 nc x ns N k=l k 1 k 


(B.9) 


Where (*) = 


Multiplying both sides of Eq. (B. 8 ) by ( q— ) 


^ *. . 9 N a * 

£2 . [ ( Q— ) a 1 + n a 1 ] = — £ cos n^, £2 ( q— ) a. 

x 52. nc ns N k=l k i £2. 


1 v *i 

7 > a k 


(B.10) 


The underlined term is the same as that in equation (B.4) 


Cancelling £2_^ on both sides, Eq. (B.10) yields 

n, . o n 0- *. 

(. q— ) a + na 1 = —• E cos nij ; 1 (. q— ) a 1 

52. nc ns N k=l k 52. k 

x x 


CB.ll) 


Equation (B.10) shows how the first time derivative term in the rotor system 

transforms into the multiblade coordinate system. It turns out that in the 

transformed multiblade coordinate system, the first time derivative term is 

to be multiplied by ( ). It can also be shown that the second time derivative 

i «i 2 

term is to be multiplied by ( ) • These multiplication factors take care 

i 

of both consistency in nondimensionalization and proper multiblade coordinate 


transformation! 

A closer look at the equation (3.11) shows that ( transforms into 

* * £2-^ i ^1 

a 1 ( 7 T- ) + na and not (a + na ) ( 7 ^— ) eventhough ( 7 ^— ) term is 
nc B, ns nc ns & 52. 

x i x 
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independent of the summation. This term does not act as a common multiplier 

for both sides. This term is multiplied only to the first derivative term 

on the left hana side. For the case when = ft , equation (B.ll) reverts 

to the original equation given in equation (A. 4) . 

Thus, in blade load expressions derived in nondimen sional form for a 

typical rotor, the time derivative terms have to be multiplied by a factor 
^1 P 

( -Q- ) (power of this factor depends on the order of differentiation) and 
i 

the same factor has to be introduced also in the time derivative terms which 
appear in the transformed multiblade coordinatv . 
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APPENDIX C 


Rotor , Blade and Body Properties 

The data provided “below describes the rotor tested in Ref. [6]. 
Rotor Geometry 

Number of blades 3 

Radius, cm 81.1 

Chord, cm 4.19 

Hinge offset, cm 8.51 

Blade airfoil NACA 23012 

Profile drag coefficient 0.0079 

Lock number 7.73 

Solidity ratio 0.0494 

Lift curve slope 2 tt 

Height of rotor hub above gimbal, cm 24.1 

Blade Mass Properties 

Blade mass (to flap flexure) , gram 209 

Blade mass centroid 

(Ref. flexure centerline), cm 18.6 

Blade flap inertia 

(Ref. flexure centerline), gram 17.3 

Blade Frequency and Damping 

Configuration 1 

Nonrotating flap frequency Hz 3.13 

Nonrotating lead-lag frequency Hz 6.70 

Damping in lead-lag (% critical) 0.52/£ 

Body Mass Properties 

2 

Rotary inertia in pitchy gram m 633 

Rotary inertia in roll, gram m^ 
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Body Frequency and Damping 


Pitch frequency, Hz 

2 

Roll frequency, Hz 

4 

Damping in roll (% critical) 

0.929% 

Damping in pitch (% critical) 

3.20% 
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